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FOREWORD 


This  is  the  final  report  on  the  Strategic  Reentry  Technology  Program 
Phase  II  (STREET-AII) .  A  two-layer  model  of  high-speed  two-  and  three- 
dimensional  turbulent  boundary  layers  with  pressure  gradient  and  surface 
mass  injection  is  presented. 

A  complete  description  of  the  turbulent  boundary  layer  computer  program, 
including  a  program  listing,  is  available  in  a  separate  AVCO/SAMSO  document. 
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ABSTRACT 


A  two-layer  model  of  the  three-dimensional  compressible  turbulent  boundary 
layer  is  developed  which  is  applicable  to  flows  with  pressure  gradient  and  sur¬ 
face  mass  transfer.  The  model  is  based  on  the  small  cross-flow  approximation  in 
which  the  spreading  metric  is  determined  by  the  inviscid  streamline  pattern.  A 
modified  Mangier  transformation  is  employed  which  permits  transformation  of  the 
boundary  layer  equations  to  a  two-dimensional  form  without  transforming  the 
turbulent  stress  or  heat  flux.  It  turns  out  that  the  computational  speed  of  the 
method  is  rapid  enough  to  enable  equations  for  the  inviscid  stream  deflection  to 
be  coupled  with  the  present  method  for  calculations  of  strong  (supercritical) 
interactions,  such  as  in  the  region  downstream  of  the  critical  point  in  re¬ 
attaching  flows  or  in  regions  of  strong  blowing. 

Solution  of  the  inner  (wall)  layer  for  the  velocity  and  enthalpy  is  ob¬ 
tained  using  mixing  length  theory  and  the  thin  layer  Couette  model.  This  solu¬ 
tion  is  obtained  without  using  a  compressibility  transformation  and  leads  to  a 
generalized  law  of  the  wall  with  mass  injection.  The  stress  and  heat  flux  in 
the  inner  layer  are  found  by  inserting  the  law  of  the  wall  into  the  full  conserva¬ 
tion  equations  and  integrating  away  from  the  surface.  In  the  outer  wake  layer, 
an  integral  moment  method  is  used  along  with  appropriate  matching  conditions 
with  the  inner  layer.  This  report  shows  that  the  dynamics  of  the  turbulent 
boundary  layer  dependsT  on  the  coupling  and  interaction  of  the  inner  and  outer 
layers.  In  fact,  the  model  is  really  designed  to  bring  out  the  nature  of  this 
inner  and  outer  layer  interaction. 

Several  solutions,  and'exp'iriraer.tal  comparisons  are  presented.  v  In  particular, 
'the  effect  of  positive  and  negative  pressure  gradients  on  the  relative  scale  of 
the  inner  and  outer  layers  is  demonstrated  as  well  as  the  effect  on  the  stagna- 
■  tion  enthalpy-velocity  relationship.  Results  are  presented  for  relaxing  flows 
where  cj  and  (2CH/Cf)  initially  are  far  from  their  equilibrium  values.  Solutions 
for  boundary  layers  far  out  of  fluid  dynamic  equilibrium  such  as  flow  over  an 
abrupt  expansion  £..id  flow  over  a  wavy  wall  are  also  given. 

The  theory  is  compared  with  a  number  of  experiments  for  both  uniform  and 
discontinuous  distributions  of  surface  mass  injection.  Solutions  have  been  ob¬ 
tained  for  values  of  M  »  pwvw/peue  as  large  as  0.05,  which  is  well  into  the 
range  of  massive  blowing.  Results  are  also  presented  for  ablation  and  trans¬ 
piration  cooling  where  the  injection  distribution  is  determined  by  an  energy 
balance  along  the  surface. 
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ay  exponent  in  eddy  viscosity  variation 

across  outer  layer,  Eq.  (46) 


Subscripts 

e  edge  of  boundary  layer 

m  at  match  point  ym 

0  isentropic  stagnation  conditions  behind  normal  shock, 

also  initial  conditions 

s  local  stagnation  conditions  at  edge  of  layer 

w  wall  conditions 

aw  adiabatic  wall  conditions 
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TASK  4.3 

A  TWO-LAYER  MODEL  OF  HIGH  SPEED  TWO-  AND  THREE-DIMENSIONAL 
TURBULENT  BOUNDARY  LAYERS  WITH  PRESSURE  GRADIENT  AND 
SURFACE  1/ASS  INJECTION. 

1.0  INTRODUCTION 


Accurate  prediction  of  turbulent  boundary  layers  about  reentry  bodies  is  of 
critical  importance  in  calculating  wall  friction  and  heat  transfer,  observables 
(initial  conditions  for  near  wake  calculations),  communication  (guidance,  and 
telemetry),  and  reentry  aerodynamics  in  general.  Recently,  interest  in  maneuver¬ 
ing  and  lifting  reentry  vehicles  has  placed  added  emphasis  on  developing  methods 
for  predicting  three-dimensional  turbulent  boundary  layers  over  bodies  at  angle 
of  attack  and  for  predicting  turbulent  boundary  layer  separation  and  reattachment. 
Separation  occurring  on  control  surfaces  and  on  the  leeward  side  of  bodies  at 
angle  of  attack  are  problems  of  much  interest,  but  are  flows  for  which  current 
theories  are  either  incomplete  or  inadequate.  For  example,  the  wake-like  model 
of  turbulent  separation  and  reattachmentl  has  yet  to  be  extended  to  three- 
dimensional  flows,  and  the  model  itself  is  incapable  of  treating  the  important 
region  downstrcm  of  the  critical  point  in  two-dimensional  reattaching  flows. 

The  objective  of  thi*  task  is  to  produce  a  method  which  will  yield  wall 
conditions  and  boundary  layer  profiles  sufficiently  accurate  for  subsequent 
calculations  of  electron  concentration  through  a  three-dimensional  turbulent 
boundary  layer.  The  work  is  aimed  at  computing  flow  fields  around  arbitrary 
bodies  at  angle  of  attack,  with  particular  emphasis  on  calculating  boundary  layers 
in  strong  positive  and  negative  pressure  gradients  with  surface  mass  injection. 
Entropy  layer  swallowing  in  regions  downstream  of  a  blunted  nose  is  also  being 
considered. 

A  two-layer  model  has  been  developed  for  this  purpose  in  which  an  analyti¬ 
cal  solution  for  the  inner  layer,  based  on  a  compressible  law  of  the  wall,  is 
matched  with  a  moment  integral  method  in  the  outer  layer.  Solution  of  the  inner 
layer  is  obtained  without  the  use  of  a  compressibility  transformation.  In  the 
outer  layer  the  explicit  appearance  of  the  density  is  removed  from  the  convection 
terms  in  the  conservation  aquations  without  transforming  the  turbulent  stress 
and  heat  flu:.. 

The  model  is  based  on  to?  tacit  assumption  that  the  inner  and  outer  layers 
can  be  treated  distinctly,  but  that  as  the  flow  proceeds  downstream  the  two 
layers  are  free  to  interact  with  one  another.  In  fact,  the  model  is  intended  to 
bring  out  the  nature  of  this  inner  and  outer  layer  interaction.  The  relative 
scale  of  the  inner  and  outer  layers  (location  of  the  match  point %,)  is  computed 
as  an  integral  property  of  the  flow,  similar  to  the  calculation  of  the  momentum 
thickness  and  displacement  thickness,  and  depends  on  the  upstream  history  of  the 
layer.  At  separation,  for  example,  which  now  can  be^predicted  by  the  model  if 
the  pressure  distribution  is  prescribed,^  .  0  ,  ym/Se  ■>  0  and  (u/uc)  -•  0  so 
that  the  boundary  layer  is  dominated  by  the  outer  layer  (or  dominated  by  the 
"wake  component"  in  the  terminology  of  Coles^)..  One  of  the  unique  features 
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of  the  model  Is  that  the  stability  of  the  system  of  ordinary  differential  equa¬ 
tions  can  be  used  to  establish  upper  and  lower  bounds  on  the  turbulent  production 

integral f rdu  in  the  outer  layer.  Stable  solutions  up  to  Mach  10  for  a  wide 
range  of  wall  to  free  stream  temperature  ratios  show  that  compressibility  has 
little  effect  on  the  properly  normalized  lateral  stress  profiles,  confirming  the 
results  of  Maise  and  McDonald. 3 

The  computational  speed  of  the  method  is  sufficiently  rapid  that  the  whole 
boundary  layer  flow  field  on  bodies  at  angle  of  attack  can  be  calculated  (except, 
of  course,  the  leeward  side  at  large  angles  of  attack).  It  also  appears  that 
the  computational  speed  is  rapid  enough  to  permit  calculations  of  supercritical 
boundary  layer-inviscid  stream  interactions. 

The  model  is  based  on  the  small  cross-flow  approximation  in  which  the 
boundary  layer  equations  are  solved  along  inviscid  streamlines.  In  this  coor¬ 
dinate  system  the  cross-flow  component  of  velocity  in  the  layer  is  uncoupled 
from  the  streamwise  momentum  and  energy  equations  and  can  be  ignored  as  long  as 
the  approximation  is  valid.  For  flow  around  bodies  at  angle  of  attack,  inviscid 
streamline  patterns  and  the  spreading  metric  are  computed  by  integrating  the 
inviscid  momentum  equation  normal  to  inviscid  streamlines  on  the  body  using 
either  experimental  or  calculated  inviscid  pressure  distributions. 4 
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2,0  TURBULENT  BOUNDARY  LAYER  MODEL 


2.1  EQUATIONS  OF  MOTION  AND  COORDINATE  SYSTEMS 


If  the  boundary  layer  equations  are  written  following  inviscid  streamlines 
and  the  cross  flow  in  the  boundary  layer  normal  to  the  local  inviscid  stream¬ 
line  direction  can  be  assumed  small,  the  mean  flow  equations  are 
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where  r(x*)  is  the  spreading  metric  determined  by  the  inviscid  flow  streamline 
geometry.  Here,  because  the  cross-flow  component  of  velocity  is  assumed  small, 
the  uncoupled  cross-flow  momentum  equation  is  ignored  and  it  is  further  assumed 
that  the  additional  contributions  to  turbulent  stress  and  heat  flux  resulting 
from  fluctuations  in  the  cross  flow  velocity  are  likewise  small. 

Thus, 
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Clearly,  for  a  body  at  large  angle  of  attack  these  stress  and  heat  flux 
models  will  break  down  in  a  region  on  the  leeward  side  because  strong  cross 
flows  with  additional  contributions  to  the  Reynolds  stress  would  have  to  be  in¬ 
cluded  along  with  interaction  with  the  outer  flow. 

The  explicit  appearance  of  the  spreading  metric  in  the  equations  of  motion 
can  be  eliminated  by  means  of  a  modified  Mangier  transformation.  Although  it 
turns  out  thatrCx  )  later  appears  in  a  constant  of  integration  in  the  compressi¬ 
ble  law  of  the  wall,  the  transformation  is  convenient  since  it  is  possible  to 
remove  r  from  the  continuity  equation  without  transforming  the  stress  or  heat 
flux.  Letting 
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where  r  and  q  are  still  defined  by  Eqs.  (4)  and  (5).  Solution  of  the  inner  layer 
is  obtained  directly  using  Eqs.  (8)  -  (10)  without  further  transformation  of  the 
coordinates  to  remove  the  explicit  appearance  of  the  density.  The  reason  for 
this  is  twofold.  First,  within  the  framework  of  the  present  model  for  the  inner 
layer  an  analytic  solution  can  be  obtained  without  introducing  a  compressibility 
transformation.  Thus,  no  useful  simplification  would  be  provided  by  such  a 
transformation.  Second,  Maise  and  McDonald3  have  shown  that  a  compressibility 
transformation-  leads  to  generally  poor  results  for  the  turbulent  stress  distri¬ 
butions  across  the  boundary  layer  even  for  constant  pressure,  adiabatic  flows. 
They  also  showed  poor  correlation  of  experimental  velocity  profiles  using  the 
transformation  for  constant  pressure  flows  with  heat  transfer.  More  recently, 
Lewis,  jit.  al.5  attempted  to  generalize  Coles'  transformation  but  they  also  found 
significant  deviation  between  experimental  velocity  profiles  and  profiles  pre¬ 
dicted  using  the  transformation,  especially  in  the  wall  layer.  They  showed  that 
the  discrepancy  between  profiles  increased  with  increasing  Mach  number  so  that 
for  flow  on  a  flat  plate  at  Mach  6,  for  example,  a  difference  of  20  percent  in 
velocity  profiles  can  be  expected.  Although  not  all  of  this  discrepancy  should 
be  attributed  to  the  transformation,  there  appears  to  be  little  justification 
for  using  it  in  the  wall  layer,  particularly  in  flows  with  heat  transfer  and 
pressure  gradient. 


In  the  outer  wake-like  layer,  however,  where  we  use  an  integral  moment 
method,  a  simple  compressibility  transformation  is  used  to  remove  the  explicit 
appearance  of  the  density  from  Eqs.  (8)  through  (10).  This  can  be  accomplished 
without  transforming  the  stress  and  heat  flux  by  letting 
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so  that  Eqs.  (8)  through  (10)  become 
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Here  r  and 

q  are  still 

the  physical  stress  and  heat 

fiux  defined  by  Eqs. 

(A)  and 

(5).  Eqs.  (13)  through  (15)  are  used  to  develop  the  integral  moment  equations 
for  the  outer  layer. 


With  the  equations  of  motion  and  proper  coordinate  systems  developed  for  the 
inner  and  outer  layers,  we  now  proceed  to  the  turbulence  model  used  in  the  two- 
layer  model. 


2.2  TURBULENCE  MODEL 


In  order  to  develop  a  scheme  for  integrating  Eqs.  (1)  through  (3),  or  their 
transformed  counterparts  Eqs.  (8)  through  (10),  it  is  necessary  to  have  some 
knowledge  either  of  the  relationships  between  r  ,  q  and  the  mean  flow  quantities 
p,  u  ,  H  (provided  such  relationships  exist)  or  have  additional  independent 
differential  equations  for  r  and  q .  Several  model  equations*^  ,7  for  the  stress 
and  heat  flux  based  on  the  latter  approach  have  been  developed  recently  and, 
hopefully,  these  will  provide  new  insight  into  the  way  the  stress  and  heat  flux 
respond  to  changes  in  boundary  conditions  as  the  flow  proceeds  downstream.  To 
date,  however,  all  these  approaches  using  differential  equations  for  r  and  q 
have  had  to  invoke  assumptions  for  closure  of  the  unknown  correlation  functions 
which  are  either  pure  hypothesis  or  whose  physical  basis  has  only  been  demonstrated 
in  incompressible  flows.  At  moderate  Mach  numbers,  however,  say  up  to  Me  -  1 
there  is  considerable  experimental  evidence  that  the  turbulent  boundary  layer  may 
be  divided  into  an  inner  wall  layer  and  an  outer  wake  layer,  as  in  incompressible 
flow.  In  the  wall  layer  the  properly  normalized  experimental  velocity  and  enthalpy 
profiles  have  been  shown  to  be  functions  of  the  local  wall  stress,  heat  flux  and 
injection  rate  and  except  for  the  streamwise  variations  of  these  quantities  are 
more  or  less  independent  of  the  upstream  history  of  the  layer. 8, 9  This  is  pre¬ 
cisely  the  result  predicted  by  nixing  length  theory  and  the  thin  layer,  Couette 
model,  which  leads  to  a  compressible  law  of  the  wall  and  a  Crocco  integral  for 
the  inner  layer.  For  flows  with  surface  mass  injection  Danberg  and  Squire  have 
shown  that  this  same  model  accurately  predicts  the  velocity  distribution  in  the 
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wall  layer.  In  fact,  the  experimental  evidence  for  the  Van  Driest^  form  of  the 
law  of  the  wall  for  M  =  0  and  the  Squire  form  of  this  law  for  M  ^  0  is  now  so 
strong  that  the  model  equations  for  the  turbulent  stress  either  "reduce"  to  this 
law  near  the  wall  or  are  matched  to  it. 

In  this  paper  a  first  approximation  for  the  velocity  and  enthalpy  profiles 
in  the  inner  layer  is  used  in  which  the  laminar  sublayer  and  ti  -i.isition  layer  is 
neglected  and  the  flow  is  assumed  to  be  fully  turbulent  to  the  vail.  Maise  and 
McDonald  have  shown  that  in  this  region  the  mixing  length  varies  linearly  away 
from  the  surface,  i.e.  (  =  K?  with  K  =  0.40,  for  Mach  numbers  up  to  about  five. 
Thus,  with  the  additional  assumption  that  Prt  =  1  in  the  first  approximation, 
the  stress  and  heat  flux  in  the  inner  layer  are 

t  -  p  K2y2  (d~/tfy)2  (16) 

q  —  p  (duV^y)  (dh/cfy)  (17) 

The  error  incurred  in  computing  integral  properties  such  as  0  and  the  wall 
stress  by  neglecting  the  sublayer  is  small,  except,  perhaps,  in  very  strong  nega¬ 
tive  pressure  gradients.  Second  approximation  profiles  are  subsequently  com¬ 
puted  which  include  the  laminar  contribution  to  the  total  stress  and  are  not 
restricted  to  pr  =  1  or  PrL  =  1. 

Even  for  flows  in  positive  or  negative  pressure  gradients  there  is  consider¬ 
able  experimental  evidence  that  the  law  of  the  wall  gives  an  accurate  represen¬ 
tation  of  the  velocity  profile  over  some  distance  ~m  away  from  the  surface.  If 
we  accept  the  generalized  law  of  the  wall  as  being  applicable  to  flows  with  sur¬ 
face  mass  injection  and  pressure  gradient,  the  remaining  questions  to  be  answered 
insofar  as  the  inner  layer  is  concerned  are  (a)  the  determination  of  the  lateral 
stress  and  heat  flux  profiles,  which  are  required  for  matching  with  the  outer 
layer,  and  (b)  calculation  of  the  distance  from  the  surface  yjn  at  which  the 
velocity  profile  departs  from  the  law  of  the  wall.  Following  Coles' 2  suggestion 
for  incompressible  flows,  the  stress  and  heat  flux  profiles  are  found  by  insert¬ 
ing  the  wall  law  into  the  full  momentum  and  energy  equations  and  integrating 
these  equations  away  from  the  surface.  In  this  way  the  stress  and  heat  flux 
in  the  inner  layrr  depend  on  the  past  history  of  the  layer  as  well  as  on  local 
properties,  and  are  unhooked  from  a  purely  local  determination  through  mixing 
length  and  the  local  velocity  gradient. 

The  determination  of  the  length  scale  y^  results  from  a  consideration  of 
the  matching  or  interaction  of  the  inner  and  outer  layers.  Thus,  attention  is 
turned  to  the  outer  layer  where  the  evaluation  of  the  stress  and  heat  flux  is 
much  more  questionable.  Bradshaw  concluded,  from  his  measurements  in  relaxing 
incompressible  boundary  layers,  that  the  stress  in  the  outer  layer  depends  on 
the  whole  history  of  the  layer  and  cannot  be  found  from  local  mean  flow  proper¬ 
ties  as  in  equilibrium  layers.  Thus,  according  to  Bradshaw  and  others,  the  use 
of  an  eddy  viscosity  or  mixing  length  theory  in  the  outer  layer  of  rapidly  ad¬ 
justing  flows  must  be  held  suspect  because  for  these  flows  there  is  no  simple 
relationship  between  the  turbulent  stress  and  mean  velocity  gradient.  If  this 
is  the  case,  then  the  relative  scale  of  the  inner  and  outer  layers,”^  ,  likewise 
cannot  be  determined  from  the  local  velocity  gradient,  as  is  usually  assumed 


in  conventional  finite  difference  methods, H»12  for  example,  but  must  be  found 
from  the  whole  history  of  the  flow.  In  the  present  theory  we  have  attempted  to 
avoid  using  a  simple  local  relation  for  determining  y*m  •  By  taking  a  higher 
moment  of  the  momentum  equation, *y  is  found  as  one  of  the  integral  properties 
of  flow. 

In  the  present  two-layer  model  it  turns  out  that  a  relation  for  the  stress 
along  the  match  point  y'm  is  required  along  with  a  relation  for  the  "production" 
integral  in  the  outer  layer 


(r/rm)  du 


The  stress  along  the  match  point  is  evaluated  using  an  eddy  viscosity  normalized 
by  the  constant  density  displacement  thickness,  as  was  suggested  by  the  calcula¬ 
tions  of  Maise  and  McDonald.  Thus  along  yl 

fm  *  Pm  'm  (W~)m  (18) 


where 


<  =  F  ue  5j*  ,  F  =  constant  (19) 

Thus,  while  the  evaluation  of  rm  is  subject  to  some  of  the  objections  of  "local¬ 
ness",  the  location  of  the  match  point  where  rm  is  evaluated  depend.'!  on  the 
upstream  history  of  the  layer. 

Because  the  turbulent  bourdary  layer  has  been  divided  into  two  more  or  less 
distinct  layers  which  are  free  to  interact  with  one  another,  the  present  model 
might  be  expected  to  reveal  certain  consequences  of  this  inner  and  outer  layer 
interaction  that  are  not  quite  so  apparent  in  other  models.  Indeed,  this  turns 
out  to  be  true  because  the  two-layer  model  provides  upper  and  lower  bounds  on 
the  magnitude  of  the  turbulent  production  integral  over  the  outer  layer.  If 
estimates  for  this  integral  are  made  which  are  too  small  or  too  large,  integra¬ 
tions  of  the  system  of  equations  downstream  develop  instabilities  where  n-0 
for  the  former  and  n  -  ~  for  the  latter.  Thus,  we  find  that  the  production  inte 
gral  over  the  outer  layer  cannot  be  assumed  negligibly  small  nor  can  it  be  evalu 
ated  using  a  constant  eddy  viscosity  across  the  outer  layer,  as  was  suggested 
by  Clauser.13  Both  assumptions  give  unstable  solutions,  although  the  Clauser 
assumption  is  stable  for  Mach  numbers  below  about  2.5. 

Thus  we  find  that  the  stability  of  the  system  of  equations,  which  reflects 
on  the  stability  of  the  interaction  between  the  inner  and  outer  layers,  provides 
a  means  of  evaluating  the  production  integral  over  the  outer  layer.  This  is  a 
most  useful  and  important  result  because  the  production  integral  contains  the 
influence  of  turbulent  intermittency  at  the  boundary  layer  edge  and  the  effects 
of  long-time  stress  "memory"  of  large  eddies  in  flows  which  are  far  from  equili¬ 
brium.  By  evaluating  the  production  integral  as  though  the  density  were  con¬ 
stant  across  the  outer  layer  and  by  using  a  cubic  variation  of  <  between  and 
u(  with  {  .Oat  the  outer  edge,  stable  solutions  have  been  obtained  up  to  Mach  10 
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for  wall  to  free  stream  stagnation  temperature  ratios  between  0.05  and  1.0. 

These  results  for  the  bounds  on  the  production  integral  and  the  assumptions  used 
to  produce  stable  solutions  verify  to  some  extent  the  results  of  Maise  and 
McDonald  who  showed  that  compressibility  has  little  effect  on  the  lateral  stress 
profiles  (when  normalized  bv  the  wall  stress)  up  to  about  Mach  5.  A  cubic 
variation  of  <  for  7m  <  7  <  is  also  consistent  with  the  <  variations  inferred 
by  Maise  and  McDonald^  and  Bradshaw.  M 


3,0  SOLUTION  OF  INNER  LAYER  AND  MATCHING  CONDITIONS 


Experiments  in  compressible  turbulent  boundary  layers  with  moderate  pressure 
gradients  and  surface  mass  injection  have  shown  that  even  at  hypersonic  speeds  there 
is  an  inner  (wall)  layer  in  which  velocity  and  temperature  gradients  normal  to 
the  surface  are  so  large  that  terms  involving  x  derivatives  in  the  conservation 
equations  are  generally  negligible.  The  lateral  extent  of  this  inner  layer  'yiri 
at  any  given  point  ~  is  determined  by  previous  history  of  the  layer  and  by 
local  conditions  such  as  magnitude  and  sign  of  the  pressure  gradient  or  strength 
of  injection.  In  large  positive  pressure  gradients  or  with  large  blowing  the 
extent  of  the  inner  layer  may  become  vanishingly  small,  while  in  negative  pres¬ 
sure  gradients  may  approach  the  edge  of  the  layer. 

With  mixing  length  theory  used  to  determine  the  stress  and  heat  flux,  and 
a  thin  layer  Coueute  model  of  the  flow  near  the  wall,  a  generalized  law  of  the 
wall  for  the  velocity  and  enthalpy  with  surface  mass  injection  is  obtained. 

These  velocity  and  enthalpy  profiles  are  then  inserted  back  into  the  fall 
continuity,  momentum  and  energy  equations,  which  are  integrated  away  from  the 
surface  to  give  the  stress,  heat  flux  and  mass  flux  at  the  outer  edge  of  the 
inner  layer.  This  operation  "unhooks"  the  stress  and  heat  flux  from  purely  local 
quantities  such  as  the  velocity  and  enthalpy  gradients.  The  velocity,  enthalpy, 
stress,  heat  flux  and  mass  flux  at  the  edge  of  the  wall  layer  can  be  obtained 
analytically,  and  these  expressions,  along  with  certain  auxiliary  relations  for 
normal  and  streamwise  derivatives  of  the  velocity  and  enthalpy,  are  used  for 
matching  with  the  integral  method  in  the  outer  layer.  One  of  the  key  features 
of  the  present  two  layer  model  is  that  the  integral  equations  in  the  outer 
layer  ultimately  determine  the  lateral  extent  of  the  wall  layer  (or  in  Coles' 
terminology  the  strength  of  the  wake  component) ;  it  is  not  specified  implicity 
in  advance  in  terms  of  the  eddy  viscosity  or  mixing  length  variation  as  in 
virtually  all  other  methods. 

With  the  possible  exception  of  boundary  layers  in  very  strong  pressure 
gradients  we  assume  that  a  wall  layer  exists  in  which  i  =  K~  and  the  velocity 
and  enthalpy  profiles  are  closely  approximated  by  a  Couetto  flow  model.  Neglect¬ 
ing  the  laminar  sublayer  (in  the  first  approximation)  and  setting  Prt  =  1,  the 
equations  for  the  stress  and  heat  flux  obtained  by  integrating  the  momentum  and 
energy  equations  away  from  the  surface  are 


r(y )  -  +  pw  ~  a  p  K2 72 


(20) 


q(y)  =qw  -*f(y)  +  p*~w  (H  -  H*)  a  p  K* 


Dividing  the  second  equation  by  the  first  and  integrating  the  result  gives  a 
Crocco  integral,  which  is  always  valid  in  the  inner  layer,  independent  of  pres¬ 
sure  gradient  or  surface  mass  injection. 

H/H„  =  1  +  B  (uVu*e)  (22) 
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or 


h^w  -  1  +  B(u/~e)-  A2(tr/tre)2 


where 


A2  =  [me/(l  +  me)  ]  ( He  /  H,, ) 

B  =  f2me/(l +me)]  (Hg/H,,)  (St/Cf) 


(23) 

(24) ' 

(25) 


Substituting  for  p  =pw  (hw/h  )  in  the  expression  for  r(*y)  ,  defining  Me  pw~w/pe~e 
and  integrating  from  the  edge  of  the  laminar  sublayer  to  any  point  'y  +  sTur/I<w  * 


with  £  a  u'/uj,  ,  gives 


f 


df 


£  ((  h/hw )  (Cf  +  Mf)] 


1/2 


r  «; 

(1  +  »e>  — 


1/2 


±  L  lL 


which  can  be  written,  after  changing  the  limits  of  integration 

r€ 

OW""  ■  i( 


v'n  f 


JQ  ((h/hw)  (Cf  +  Mf)l 
where  B  is  the  empirical  "constant” 


1/2 


—  fc,  7+  *  B 
K 


B  a 


-  Y  ^Ts+  +  (hw/he)~1/2  j 


d( 


1/2 


0  [(h/h,,)  (Cf  +  M£)l  i 


Whereas  in  incompressible  flow  Bfi:5,  the  experiments  of  Danberg  and  Squire  at  I 
supersonic  speeds  have  demonstrated  that  B  is  a  function  of  M  and  H^H,,  ,  which 
can  be  approximated  by  the  curve  fits*:  J 


with 


0  <  Me  <  3.5  :  B/BM_0  =  1  -  0.6  x  103  M  (Me/3.5)3/i 
Me  >  3.5  :  B/BMb0  -  1  -  0.6  x  103  M 

Mc> 5:  BMb0  -  5  +  B  (1-Hw/He>, 

H_,\  ( 


„  t 

B">0 


Me<5:  BM  =  0 


5+8  1  - 


He  )\  5 


linnet  ctVPpolatioa  of  5  data  presented  in  Danberg  was  used  for  Jl^, /If.  <0.5.  P/ofiie  data  is  desperately 
vcik/  Jcr  -rodify)  the  values  of  D  (and  values  of  a  ,  Section  6.0)  that  are  indicated  for  highly  cooled  walls. 


needed 
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I  ! 


I 


1 


Defining* 


H 

Cj  m  (B  -  VB2  +  4  A2  )  /2Kl 

! 

* 

ad 

s 

1 

i 

.  c?  s.  Cf/M,  -’k2  K  (C2  -  Cj  )i 

a  -  ’ 

1  i  ‘ 

H  • 

1  and  since  k  ?+  «=  L  y+  +  In  1,/f 

..-.2 


the  'Compressible  law  of  the  wall  becomes 

1/2 


-2 


he/hw 


A2M(C3  +  c2) 


[y(d,k)-.y(Ofs0,k)]  =  -i.  In  t(Cf')1/2  Rey  ] 


I  I 


1  .  L 

+  -  in  -  +  B 


(26) 


’where 


J  (fl’.k)  a 


f- 

Jn  (1-i 


d  6  ■ 


-q  (1  -  k2  sin2  0)  ! 


I 


y+  h  y  ur/i/w  =  .(Cf')1/2  Rey 


l1  Fn 


He_ 

Hw 


V(W) 


(3y-l) 

(l+oe!)  2<*  U  Mc  y 


(27) 


(28) 


(29) 


For  M  nO,  k  =>0and  the  elliptic  integrals  J(0,k)  reduce  to  arcsin  functions,  and 
the  law  of  the  wal}.  is  glyen  by  the  Van  Driest  "generalized"  velocity 
distribution. 10  Evaluating  the  law  of  the  wall  at  ym  gives  the  velocity  ratio 
(u / “e )  a  Um  along  the  match  point.  Actually,  what  we  require  for  matching  with 
the  outer  layer  is  the  streamwise  derivative  .of  (u/ue)  ,  which  is  found  by 
•  differentiating  the  law  of  the  wall.  f1 

'  i  1  ' 

With  the  velocity  profile  given  by  the  law  of  the  wall  and  the  enthalpy 

profile  given  by  the  Crocco  integral  in  the  inner  layer,  exact  expressions  for 

the  stress  and  heat  flux  at  ym  can  be  obtained  by  integrating  the  momentum 

and  energy  equations  away  from,  the  surface,  i.e.,, 


•For  Urge  blowing,  sin2  0  and  k2  are  redefined,  i.e.  k2  ■  (C}  +  C2)/(C2  -  Cj)  and  gin2  0  m  (C2  -  f)/(C2  +  Cj)  no 
that  k2  is  always  less  than  one. 
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I 


I 


where  the  integrals  in  these  expressions  are  evaluated  using  Eqs.  (23)  and  (26). 
For  example,  if 


which  can  be  integrated  repeatedly  by  parts  if 
a  «  K(ue/ur),  Z  «»  csp  [ct  f(  £)  ] 
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where  — '*-0(Cf1/2)  <  <  1  .  By  approximating  h/hw  -  hm/hw  in  f(£)  and  g(f)  and 


1  «  ^  *  l,JlUm  in  f(£)  the  above  series  truncates  with  the  term  0(a-1) 


for  Jj  ,  and0(a-2)  for  J2  so  that* 


K  l  u. 


W  <VM* 


M  - 
1  +  —  Un 


The  expressions  for  Jj  and  J2  are  the  key  results  for  the  evaluation  of  the 
stress,  heat  flux  and  mass  flux  in  the  inner  layer.  With  Jj  and  jy  given  by 
Eqs.  (33)  and  (34)  the  evaluation  of  rm  and  (q  +  u  r  )m  is  straightforward 
(although  somewhat  lengthy)  and  the  results  are  given  in  Appendix  1. 


Solution  of  wall  layer  for  the  velocity,  enthalpy,  stress  and  heat  flux  has 
been  obtained  analytically  in  terms  of  the  given  wall  and  edge  conditions 
(Me,Hw/He  and  M  );  three  "parameters"  (ym,  C(  and  St  ),  and  two  empirical 
"constants"  (K  and  B).  By  matching  this  inner  solution  with  an  integral  moment 
method  for  the  outer  layer  the  three  "parameters"  ym  ,  C[  and  St  are  determined. 
In  addition  to  u/oe ,  H/He  ,  r  and  q,  du/dy  and  dH/dy  are  also  matched.  These 
expressions  are  given  in  Appendix  1  along  with  certain  other  auxiliary  expres¬ 
sions  (such  as  d£Tm/dx,  dHm/dx,  etc.)  required  for  the  matching. 


*  These  approximations  yield  identical  results  for  the  leading  term  of  the  scries  as  the  exact  expressions  for  f  and  g, 
and  are  used  here  merely  to  simplify  the  task  of  differentiating  J2 . 
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4.0  MOMENT  INTEGRAL  METHOD  FOR  OUTER  LAYER 


In  order  to  determine  the  relative  scale  of  the  inner  and  outer  layers,  ym , 
and  the  quantities  rm  and  qm  ,  which  are  related  to  the  quantities  of  interest, 
namely  C(  and  St ,  in  the  equations  given  in  Appendix  1,  a  moment  integral  method 
is  used  in  the  outer  layer.  In  this  layer  independent  two-parameter  families  of 
profiles  are  chosen  for  u  and  H.  The  "free"  parameters  of  the  profiles  are 
determined  implicitly  by  matching  with  relations  obtained  for  the  inner  layer 
(law  of  the  wall  and  Crocco  integral)  for  u  -and  H  and  .'he  lateral  and  streamwise 
derivatives  of  these  quantities  along  the  unknown  match  point  ym  .  The  profiles 
selected  are  power-laws  of  the  form 

U  a  U/Ue  =  1  -  (1  -  Um)  (1  -  T))n  ,  n  >  0 
H  a  H/He  =  I-U-HjjXI-ijJ®,  m  >0 

so  that  very  "full"  profiles  without  overshoot  are  possible  with  n,  m>>l  .  In 
order  for  a  Crocco  integral  to  exist  in  the  outer  layer  m  must  equal  n  but  in 
general  this  will  not  be  true.  In  fact,  in  terms  of  the  more  conventional 
definitions 


(35) 

(36) 


► 


CH  ■  -qw/(Hw-He)  pe  ue 


2  me  St 

(1  +me)(l  -  H^/Hg) 


(37) 


~f  a  2t w/pe  Ue2  =  2  Cf 


(38) 


it  is  easy  to  show,  using  Eq.  (A  1-2)  and  the  definition  of  B  that  a  Crocco 
integral  (m  ■= n)  exists  in  the  outer  layer  only  when  the  Reynolds  analogy  holds, 
i.e.,  when  CH  =  q/2  . 


In  the  above  profile  families 
V  ■  <V-Ym>/« 


(39) 


where  8  is  the  thickness  of  the  outer  layer  in  the  transformed  coordinate 


Y  - 


ae  P 


aiPo 


dy 


(40) 


Because  the  momentum  and  energy  equations  are  satisified  exactly  in  the  inner 
layer  using  the  law  of  the  wall  (Eqs.  (Al-5)  and  Al-6)  ),  no  new  information 
would  be  obtained  by  integrating  the  conservation  equations  and  velu-lty  moment 
across  the  inner  layer.  The  momentum  integral,  first  velocity  moment  and  energy 
integral  for  the  outer  layer  are  obtained  by  integrating  the  transformed  con¬ 
servation  equations  from  Y  =  Ym  to  Y  =>  Ym  +  8 . 
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Momentum  Integral  -  Outer  Laver 


dU, 


*  .  «■  .  dn  dS 

1  dx  ■  +  SA2  17  +  H  +  (P,/P0H l+me) 


-(y+  D 
2 (y-i) 


•  M  (i-  um) 


(41) 


1  dr 


—  .—  I  .  (E,-!,-,,,-!  -1^  .  / 

\PoUe  dx  /  1  5  2'  M  d 

'  *  \V*ePo 

First  Velocity  Moment  -  Outer  Laver 


m 


d  U_ 


—  (y+1) 


~«a3  dx  +  SA4  —  +  I6  ii  +  (ps/p0)(l  +me)  2(y_1)  •  M(l-Um2)J 


-  (1- 


dm 

dx 


5  dMc 
-16-2,3)—  — 


(42) 


Energy  Integral  -  Outer  Laver 


-2 


V  *e  Po 


|  rm  <Um  +  <\>> 


d  U  . 

e„  m  _  dn 

-SB*  - -  +  SB 


dS 


dHm 


dm 


dx 


—  +  E4  +  SBj  —  +  SB,  — 

dx  dx  *  dx  1  dx 


-(y+D 

+  (ps/p0)(i+me)  2(y_1)  •  M  (i-I^)  -  (i-  TL)  f— —  ~ 

1  P0  Ue  dx 


8  dMe 

~  E<1  Me  dx 


/o  ueHey 


q“  ~  *  P  a  H  1  Um  r“> 

ro  O  nC 


(43) 


Jnte8^1  functions  A;,  B;,  Ej  and  I;  are  given  in  Appendix  2.  The 
equations  for  m/p0Ue  and  the  normalized  stress  and  heat  flux  at  ym  ,  obtained 

t  JInr  Jay!r>  are  8iVSn  in  APPendix  The  term  involving  the  integral 
of  the  turbulent  stress  across  the  outer  layer  in  the  velocity  moment  (the  & 
turbulent  production)  is 
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SSSOSonsnacwe* 


fi 


(44) 


According  to  the  work  of  Maise  and  McDonald,  there  is  virtually  no  effect  of 
compressibility  on  the  normalized  otress  distribution,  r/rm  ,  across  the  turbu¬ 
lent  boundary  layer.  Thus,  we  can  evaluate  r/rm  as  though  the  density  is  con¬ 
stant,  i.e., 


r  *  /  e  \  (  dV 
/  (dU 


(dU/d.?) 


(dU/d^)m 


(45) 


Since  t  0  at  the  edge  of  the  layer,  the  lateral  variation  of  f  in  the  outer 
layer  can  be  approximated  by 


1  -«?" 


(46) 


where,  according  to  measurements  by  Bradshaw  and  the  analysis  of  data  for  super¬ 
sonic  layers  by  Maise  and  McDonald,  o>  has  a  value  around  2  or  3.  On  the  basis  of 
stability  of  solutions  over  a  wide  range  of  Mach  numbers  we  tentatively  select 
the  value  m  =  3 .  (See  discussion  in  Sections  2  and  8.) 

Finally,  a  relation  for  r  is  obtained  in  terms  of  the  eddy  viscosity  and 
velocity  gradient  at  the  matching  point  ym .  Following  the  suggestion  of  Herring 
and  Mellor  and  the  results  of  Maise  and  McDonald,  the  length  scale  for  the  eddy 
viscosity  is  the  constant  density  displacement  thickness 


8:* 


f 

»  A 


H  -  u/u..)  dy' 


(47) 


and 


<m  -  F  ue  Si  •  F  ~  °-018- 
Then,  since 


(48) 


'm 


Pm  fm 


and(duVdy')  if  known  from  the  solution  of  the  inner  layer,  the  normalized 


stress  along  ym  is  given  by 
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_  (y+i) 

=  (p*/p0>2  +  ®e>  Y  1 


F  •  (Sj*/5)  a*  n 
tHm  +  »e  (Hn,  ~  Um2  )  ]  2 


(49) 


With  the  three  integral  equations  for  the  outer  layer,  the  above  equation 
for  rm  ,  and  the  matching  relations  from  the  inner  layer,  the  system  of  ordinary 
differential  equations  can  be  integrated  downstream.  Since  the  solution  of  the 
inner  layer  involves  three  "parameters"  (Cf>  St  and  ym),  and  the  unknown  in  the 
outer  layer  is  5 ,  these  four  quantities  must  be  specified  at  the  initial  station 
in  order  to  form  all  other  quantities  (n,  m,  etc)  <and  all  the  initial  streamwise 
derivatives. 
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5.0  INITIAL  CONDITIONS  AND  INTEGRAL  THICKNESS 


Instead  of  having  to  specify  ym  and  8  along  with  Cf  and  St  at  the  initial 
station  *=0,  it  is  more  convenient  to  specify  0  and  8e  ,  the  momentum  thickness 
and  thickness  of  the  layer.  From  the  equations  for  Jf  and  J2  in  the  inner  layer 
and  the  integral  functions  in  the  outer  layer  one  obtains  the  following: 


y  +  1 

St  “  ym  +  (po/p.W+ne)  2(V  15  Id +«e)Ei  -  me  I2  ]  8 


a  UB0-Un)ym  C{1'2  ym 

0  »  — - - -  + 


(hm/he>  K(hm/be)1/2 

y  +  i 

+  +me)  1  Eli  —  I2  1  8 


(2Um-l)/30-(2/K2)Cfy  j8( 


m 


(50) 


(51) 


8* 


U„ 


1  - 


(h«/he) 


Cf1/2  yB  Po 

yB  +  - ttt 

K(hm/he)l/2 

y  +  1 


+  (P0/P* )  d  +n>e)  2(y_1)  ((1  +  n>e)E1  -  I2  —  It  ]8 


V  °  yB(1-um)  +  (Cfha/he)1/2  (ym  /30/K) 

t  y+1 

/ Po\  2 (y  -  1)  r 

+  1—1  d+  me)  r  [  (1  *  me)  (Ej  -  E2)  -  me  (I2  —  Ij)  ]  8 


Also, 


(52) 


(53) 


8t  -  (r/L)  8*  ■  (t/L)  8* ,  0-(r/L)7 


(54) 


By  specifying  Cf  ,  Cjj,  8e  and  <1  at  the  initial  station,  the  above 

matching  conditions  are  used  to  calculate  the  remaining  unknowns: 

n  ,  H  ,  and  B  . 
m 


expressions  and 
ym  » um  »  8  » m  » 
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6.0  FINAL  VELOCITY  AND  ENTHALPY  PROFILES 


Since  at  hypersonic  speeds  the  peak  boundary  layer  temperature  can  occur 
well  within  the  wall  layer,  accurate  velocity  and  enthalpy  profiles  near  the 
wall  are  required.  In  order  to  avoid  the  logarithmic  singularity  at  the  surface 
given  by  the  law  of  "the  wall,  the  laminar  sublayer  and  transition  layer  are  in¬ 
cluded  in  the  determination  of  final  profiles  (as  well  as  arbitrary  laminar  and 
turbulent  Prandtl  numbers)  so  that  u'=  0  at  *=  0  .  This  is  accomplished  by  in¬ 
cluding  the  laminar  contribution  to  the  total  stress  and  heat  flux  in  the 
expressions  for  r(y')  and  q  (y') in  the  wall  layer  and  by  using  the  Van  Driest 
variation  for  the  dimensionless  mixing  length  <t>  ,  i.e., 

♦  =  K'y+  ( 1  -exp(-7+/a>]  (55) 

where 

t(7f .  «)  =  ~r  1/ * w  •  (56) 

For  adiabatic  flows  with  M  =  0 ,  a  is  a  constant  and  equals  the  incompressible 
value  a  n  25  .  For  cold  walls_and  M  ^  0  the  experimental  velocity  profiles  of 
Danberg  and  Squire  show  that  B  decreases  with  increasing  injection  and  indicate 
large  variations  in  a  .  The  determination  of  a  from  these  experimental  profiles 
is  discussed  below. 

Introducing  the  normalized  wall  layer  variables 
**  -  ~/~r-  h+2  h\,  7+  *  ~r  y/vw 

the  differential  equations  for *u+  and  h+  obtained  from  the  equations  for  the 
total  stress  and  heat  flux  become 


with  boundary  conditions,  ~+ =  0  : ~+  *  0  ;  h+  =  1  . 

The  differential  equations  for  *u  +  and  h  +  need  only  be  solved  at  selected 
points  in  the  flow  field  where  accurate  profiles  near  the  wall,  including  the 
sublayer,  are  required.  Since  both  boundary  conditions  for  the  two  first  order 
equations  are  known  at  the  wall,  the  solution  does  not  require  iteration  and  is 
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performed  as  an  auxiliary  calculation  in  the  downstream  integration  of  the  two- 
layer  equations. 

The  asymptotic  integral  of  the  equation  for'u+  for~+>  >  1  is 


_ _ 

[h+(l+vw-tr+/rr)l1/2 


(59) 


which,  for  a  given  value  of  a  and  solution  of  the  differential  equations  for  u’+ 
andh+,  can  be  solved  for  B.  By  iterating  on  an  assumed  value  of  a  and  matching 
the  computed  B  with  the  experimental  values  reported  by  Danberg  and  Squire  the 
variation  of  a  with  Hw/He  and  M  for  various  Mach  numbers  was  obtained.  (See 
Figure  1.) 


Figure  1  a  VARIATION  WITH  Me,Hw,vw  OBTAINED  FROM 
EXPERIMENTAL  VALUES  OFl 


-22- 


Eqs.  (57)  and  (58)  determine  *u  +  and  h+  as  functions  of  y’+  in  the  inner 
layer,  i.e.,  0  <  y'+  <  yj  .  The  profiles  of  'u/u'e,  h/he  and  H/He  as  functions 
of  ~  are  found  from 


u/7e  =  *u  +  [(hw/he) Cf ] 
h/he  =  h+  (hw/he) 

H/He  =  (1  +  mc)-1  [(h/he)  +  me(u'/ue)2] 


(60) 

(61) 

(62) 


where 


-1 


y  =  y 


11 

So  J 


Po  M 


(3y- 1) 


Me(l  +  me) 


^  He  V  +  a 

d  ,c'', 


1/2 


(63) 


In  the  outer  layer  u’/u'e  *  U  and  H/He  «  H  are  calculated  using  Eqs.  (35)  and 
(36)  and  (h/he)  is  computed  from  Eq.  (62).  The  distance *y  is  found  from  Eq.  (11). 


=  Vrr 


o  (t)s  j 


y 


+ 


(h/hc)  d  7/ 


(64) 


7.0  COMPUTER  RESULTS  AND  EXPERIMENTAL  COMPARISONS 


7.1 

By  specifying  initial  values  for  5e  ,  0 ,  and  CH,  and  also  the  unit 
stagnation  Reynolds  number,  where  if  7.° 

I  o-r±- 

P0*Jp0  =  (Ue/i'e)  I  (Ps/Po^  Mr  (1  +  me)  ^  ^ 

Eqs.  (23),  (26),  (50),  (51),  (a  1-1)  and  A  1-  2)  are  solved  simultaneously  to 
give  all  the  necessary  starting  conditions.  A  set  of  nine  nonlinear  differential 
equations  is  then  integrated  to  give  the  solution  downstream  of  the  initial 
station.  These  differential  equations  are  the  momentum,  first  moment  and  energy 
equations  for  the  outer  layer,  the  equation  for  rm  from  jthe_inner  layer,  and 
equations  obtained  by  differentiating  the  relations  for  Um,  Hm>  n  and  B  with 
respect  to  x.  The  method  used  to  solve  the  system,  which  is  linear  when  solving 
for  the  derivatives,  uses  Gauss  elimination  with  partial  pivoting,  with  an 
iterative  routine  to  improve  the  solution. 15  a  comparison  of  this  method  of 
solving  systems  of  linear  equations  with  other  techniques  is  given  in  Reference 
16. 

For  zero  pressure  gradient,  nonadiabatic  flows  the  average  computing  time 
on  an  IBM  360-75  machine  is  about  15  seconds  for  Re#  to  increase  by  a  factor 
of  102.  This  is  approximately  the  computing  time  required  by  finite  differ¬ 
ence  methods  to  calculate  a  single  station.il  As  of  the  date  of  completion 
of  this  report,  computer  rur.a  have  been  made  for  adiabatic  and  cold  walls 
(0.05  <  <  1)  up  to  Mach  10,  relaxing  flows  where  Cf  and/or  CH  are  far 

from  their  equilibrium  values  initially,  boundary  layers  in  strong  positive  and 
negative  pressure  gradients  and  flows  with  a  varying  spreading  metric.  Compari¬ 
sons  have  been  made  with  several  of  the  experiments  for  incompressible  flows 
reported  at  the  AFOSR-Stanford  conference^-2  and  results  for  £ff ,  0  and  S*/0  com¬ 
pare  quite  well  with  most  of  these  experiments.  For  example,  Figure  2  shows  a 
comparison  with  one  of  the  more  interesting  experiments  —  the  relaxing  flow 
from  a  positive  pressure  gradient  reported  by  Bradshaw  and  Ferriss.  Momentum 
balance  considerations  show  that  the  first  experimental  value  for  (at x  B  4.4  ft.) 
is  probably  too  small.  If  the  initial  0  had  been  taken  to  be  about  0.475  in.  the 
theoretical  curve  would  have  just  about  passed  through  the  five  values  of  6 
measured  farther  downstream. 

Figure  3  shows  that  Eqs.  (57)  and  (58)  yield  quite  accurate  results  for  the 
velocity  profile  in  the  wall  layer  with  surface  mass  injection,  provided  that 
emprical  values  of  the  additive  constant  d  are  available.  These  results  show 
that  not  only  is  the  linear  law  of  the  wall,  (Eq.  (59)),  still  in  evidence  with 
blowing  but  that  the  value  of  y+  where  the  profile  deviates  from  a  linear  law 
can  be  represented  reasonably  well.  Here  f(u+)  is  the  integral  on  the  left  hand 
side  of  Eq.  (59). 

The  variation  of  the  equilibrium  flat  plate  cf  with  Reg  and  Me  is  compared 
with  several  experiments  in  Figures  4  and  5.  These  results  (with  the  exception 
of  the  experiments  by  Lee,  et. _al. )  are  for  adiabatic  flow,  however,  the  theory 
predicts  a  relatively  weak  effect  of  wall  cooling  on  skin  friction.  At  Mach  6 
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and  Reg  =10^,  for  example,  Cj  for  !HW/He  =  0.05  is  only  about  15  percent  greater 
than  the  adiabatic  value  of  Cf .  At  these  same  conditions  the  Reynolds  number 
based  on  8*  for  a  cold  wall  is  about  one-third  the  value  for  an  adiabatic  wall. 

I  :  . 

Figure  6  demonstrates  the  existence  of  an  equilibrium  Cf  variation  with  Reg 
and  the  approach  to!  this  equilibrium  variation  if  the  initial  cTf  is  too  large  or 
too  small.  A  similar  typd  of  relaxation  for  the  heat  transfer  and  total  enthalpy 
profile  is  shown  in  Figure  7.  For  these  solutions  the  initial  Cf  was  set  equal 
ito  tjie  equilibrium  value  but  the  initial  value  tof  the  Reynolds  analogy  parameter 
■  2  CH/£f  was  set  at  values  other  than  oive.  This  resulted  in  an  initial  H  versus 
u  profile  which  Was  nonlinear  in  the  outer  layer.  The  solutions  shown  in  Figure 
7  demonsttate  the  ^ela^ation  of  the  Reynolds  analogy  parameter  back  to  one,  and 
relaxation  of  an  initially  perturbed  H  versus  u  profile  back  to  a  linear  Crocco 
integral  for  the  whole  layer.  (See  schematic  inserts  of  profiles  in  Figure  7.) 
Relaxing  flows;  of  this  ;sort  have  been  measured  on  wind  tunnel  walls  downstream 
of  the  nozzle  throat.  i 

i 

Figurd  8  shows  results  for  Cf ,  6t  Um,  and  the  Reynolds  analogy  parameter 
■for  a  boundary  layer  in  a  region  where  dp/dx> >  0  followed  by  a  region  where 
dp/dx  <  0  l.  The  edge  Mach  number  decreases  from  2.6  to  1.3  over  the  first  50 
inches  and  then  increases  to  2.6  again  atx  =  100  inches.  In  each  case  the 
results  are  compared  with  a  solution  for  Me  =  constant  “2.6,  In  the  region  of 
positive  pressure  gradient  the  wall  layer  thickness  decreases  but  then  increases 
as  dp/dx  <  0.  This  behavior  is  reflected  in  the  variation  of  um,  which  shows 
that  hear  x  =  100  inches  most  of  the  velocity  variation  across  the  boundary  layer 
occurs  in  the  wall  layer.  Thus,  the  two- layer  model  is  seen  to  predict  the  kind 
of  behavior  hypothesized  by  Coles^  in  his  paper  on  the  law  of  the  wake. 

I  !  • 

Figures  9  and  10  show  comparisons  of  the  theory  with  datal8,19  f0r  two-dimen¬ 
sional  flows  in  negative  pressure  gradients,  at  initial  Mach  numbers  of  1.5  and  3.9, 
respectively,  while  Figure  11 -shows  a  comparison  with  the  method  of  Bradshaw  for 
a  Iflow  in  zeto  pressure  gradient  at  Mach  3,  followed  by  a  region  of  strong 
adverse  pressure  gradient.  Both  methods  predict  that  the’flow  separates,  i.e., 

Cf  -*0.  For  the  two-layer  model  the  solution  was  stopped  when  Cf  decreased  to 
2  x  10“5  at:  x  =■  67.8  cm.  In  this  region  Um  ■•♦  0  and  the  layer  was  composed 
almost ( entirely  of  the  outer  layer,  a^ain  demonstrating  the  kind  of  behavior 
predicted  by  Coles. 

Comparison  with  data  obtained  by  Lewis'^  for  boundary  layers  in  an  initially 
positive  followed  by  a  negative  pressure  gradient  are  shown  in  Figures  12  and  13. 

In  the  "strong"  pressure  gradient  case,  the  Mach  number  varied  from  A  to  2.57  to 
A, 17  over  a  distance  of  about  3A  inches.  In  the  "moderate"  pressure  gradient 
case  the  Mach  number  varied  from  A  'to  3.2  to  3.8.  The  experiments  were  performed 
with  tn  outer  hollow  cylinder  and  an  inner  axisymmetric  pressure  generating  body. 
The  boundary  layer  data  was  taken  ion  the  inside  surface  of  the  outer  cylinder  in 
brder  to  eliminate  effects  produced  by  longitudinal  curvature.  The  curves  for 
S*  and  0  are  compared  with  predicted  variations  for  zero  pressure  gradient  in 
both  Figures  12  and  13.  The  ],ocal  decrease  in  ?*  in  the  region  of  strong  posi¬ 
tive  pressure  gradient  (Figure  12)' demonstrates  the  supercritical  behavior  of 
turbulent  boundary  layers  at  high  Mach  numbers.  Figure  1A  shows  the  predicted 
skin  friction  variations  for  the  two  cases.  Despite  the  relatively  large  peak 
pressure  atlo  of  about  8  in  the  strong  pressure  gradient  case,  the  theory  pre¬ 
dicts  th.  t  the  skin  frictidn  coefficient  is  reduced  by  only  about  50  percent 
from  the  local  flat  plate  Cf.  , 
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Figure  6  RELAXING  FLOWS  AT  Me  =  CONSTANT  =  2.6 


Me  ■  3.6,  Hw/He  »  0.5 
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1 1  COMPARISON  OF  TWO-LAYER  MODEL  WITH 
BRADSHAW  METHOD,  Hw/He  =  0.9,  Re0Q  =  103 
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Figure  13  COMPARISON  WITH  LEWIS,  MODERATE  PRESSURE  GRADIENT 


Figure  15  shows  several  solutions  (various  initial  values  of  cf  )  for  the 
experiments  conducted  by  McLafferty  and  Barber, 21  which  were  performed  on  a  two- 
dimensional  curved  compression  ramp.  The  Mach  number  decreased  from  3  to  1.8  in 
a  distance  of  3  inches,  and  the  boundary  layer  was  tripped  upstream  of  the  ramp 
by  normal  injection.  Thus,  initially  was  probably  slightly  less  than  the 
equilibrium  value  of  1.8  x  10“3.  Figure  15  shows  that  for  an  initial  Cjp  ■  1.75  x 
10-3  the  flow  remains  attached  (McLafferty  and  Barber  reported  that  the  flow  did 
not  separate)  but  for  an  initial  Cf  =  1.5  x  10”  3  the  flow  separates  at  x  »  2.73 
inches.  In  the  latter  case  the  strong  pressure  gradient  did  not  permit  relaxa¬ 
tion  back  to  the  equilibrium  Cf  so  that  flow  separation  on  the  ramp  could  have 
been  produced  if  the  injection  upstream  of  the  ramp  had  been  stronger. 

Comparison  of  the  theory  with  data  taker,  on  an  axisymmetric  curved  compres¬ 
sion  ramp  is  shown  in  Figures  16  and  17.  In  this  case  the  Mach  number  decreased 
from  5.75  to  2^6  and  the  spreading  metric  increased  by  a  factor  of  1.7.  Figure 
17  shows  that  Cf  actually  increases  slightly  due  to  the  metric  effect  (the  two- 
dimensional  results  show  only  a  40  percent  decrease  in  Cf  at  the  trailing  edge) 
and  the  layer  does  not  even  come  close  to  separating.  Both  theory  and  experiment 
show  a  substantial  increase  in  the  wall  heating  (Figure  17). 

Comparisons  with  data  on  a  waisted  body  of  revolution  at  M*,  =  1.4  and  2.4 
are  shown  in  Figures  18  and  19  which  indicate  large  effects  produced  by  the 
spreading  metric,  particularly  on  the  momentum  thickness  variations. 

Figure  20  shows  some  results  for  relaxing  flows  produced  by  taking  a  given 
initial  turbulent  boundary  layer  profile  and  letting  it  relax  downstream  over 
conical  bodies  with  different  cone  angles.  These  solutions  illustrate,  for 
example,  the  behavior  downstream  of  a  sphere-cone  junction.  The  solution  for 
/3  »  90  degrees  is  representative  of  the  flow  over  the  windward  side  of  a 
sphere-cone  at  large  angle  of  attack,  i.e.,  the  effective  cone  angle  /3  a  0C+ a 

is  90  degrees,  where  0C  is  the  cone  half  angle.  The  solutions  were  all  obtained 

for  a  boundary  layej^  initially  in  equilibrium  for  a  constant  metric  (flat  plate 
or  cylinder),  with  80/r0  =  0.1.  Thus,  the  solution  for  /3  =  0  degrees  remains 

in  equilibrium  for  7  >  70  .  The  results  for  /3  >  0  degrees  all  show  that  the 

momentum  thickness  grows  more  slowly  with  7  than  if  /3  =  0,  while  the  skin  fric¬ 
tion  coefficient  increases  with  increasing  /3  at  a  given  point  downstream  of  70. 
In  fact,  for  large  /9  (see  results  for  (i  =  90  degrees)  the  momentum  thickness 
actually  decreases  at  the  beginning  of  the  relaxation.  Far  downstream  the 
curves  for  Rcq  (and  also  Cf)  for  /3  >  0  degrees  all  coalesce  into  a  single 

curve,  which  is  displaced  from  the  ft  =  0°  solution  by  a  constant  factor  in 
Re(7-7 ) •  Table  I  shows  results  obtained  for  the  various  values  of  /3  at 

Re0  ■  3  x  10^  from  which  it  is  evident  that  the  boundary  layers  for  /3  >  0° 

are  virtually  identical  to  the  boundary  layer  for  fi  “  0°  if  Re/~  ~  \  for 

/3  >  0°  is  about  twice  the  value  of  ~  for  p  *»  0°.  This  is  a  well  known 

result  for  flow  over  sharp  cones  and  flat  plates  in  equilibrium. 24 

Some  results  computed  for  a  turbulent  boundary  layer  around  a  sphere-cone 
nosetip  at  a  «=  0  and  20°  angle  of  attack  (windward  side)  are  shown  in  Fig¬ 
ures  21a,  b,  c.  The  solutions  were  performed  by  starting  15  degrees  from  the 
stagnation  point  with  a  transition  Rcq  ■  200  and  several  values  of  . 
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COMPARISON  WITH  DATA  FOR  dP/dx  >0  Figure  16  COMPARISON  WITH  DATA  ON  AN  AXISYMMETRIC 

COMPRESSION  SURFACE 


For  these  conditions  the  solutions  became  independent  of  the  starting  value  of 
C(  at  about  35  degrees  from  the  stagnation  point.  The  solutions  demonstrate  the 
relatively  slow  growth  of  6  and  Se  along  the  windward  side  at  angle  of  attack. 
Along  the  windward  meridian  increases  by  about  35  percent  as  a  increases 
from  0  to  20  degrees.  Also,  2C(j/Cf  is  less  than  one  in  both  cases  which  is 
characteristic  of  turbulent  boundary  layers  in  negative  pressure  gradients. 

|  (See  Figures  7  and  8c.) 

|  TABLE  I  ' 


EFFECT  OF  CONE  ANGLE  FOR  A  GIVEN  INITIAL  PROFILE 
(SOME  RESULT  SAT  EQUILIBRIUM  FAR  DOWNSTREAM  OF  INITIAL 
STATION.  Me  =  2.6) 


r 

Re/'*'  **'  . 
(x-x0) 

cf 

Re? 

Res'. 

0 

3.81  x  107 

1.314  x  10~3 

3.02  x  104 

1.22  x  105 

15 

8.35  x  107 

1.323  x  10-3 

3.04  x  104 

1.22  x  105 

30 

8.49  x  107 

1.322  x  10-3 

3.06  x  104 

1.23  x  105 

60 

8.41  x  107 

1.325  x  10-3 

3.02  x  104 

1.22  x  105 

90 

8.45  x  107 

1.324  x  10-3 

3.03  x  104 

1.22  x  105 

In  Equilibrium,  Boundary  Layer  on  a  Cone  is  Equivalent  to  that  on  a 
Plate  if  Cone  Rej.~  2.2  x  Plate  Rex 

Figure  22  shows  some  results  for  as  a  function  of  Re^  calculated  for 

relaxing  flows  at  a  constant  Mach  number  of  9.37.  Several  solutions  with  various 
starting  values  of  c’f  and  2  Ch/Cj  were  performed  to  obtain  a  qualitative  com¬ 
parison  with  experiments  performed  by  Laderman  and  Demetriades^  in  the  boun¬ 
dary  layer  on  the  wall  of  a  wind  tunnel.  The  various  initial  conditions  were 
used  to  try  and  ascertain  the  effect  of  the  strong  negative  pressure  gradient 
at  the  nozzle  throat  on  the  downstream  properties  of  the  boundary  layer.  The 
experimental  profile  data  of  Laderman  and  Demetriades  was  taken  at  a  station 
160  inches  downstream  of  the  nozzle  throat  and  shows  a  nonlinear  stagnation 
enthalpy  versus  velocity  profile  characteristic  of  turbulent  boundary  layers  on 
wind  tunnel  walls.  This  nonlinear  profile  has  been  interpreted  as  evidence  that 
the  layer  is  still  in  a  state  of  relaxation  produced  by  the  strong  negative 
pressure  gradient  at  the  nozzle  throat. 

The  theoretical  solutions  were  started  at  an  arbitrary  value  of  Re^,  °  1,5  x 

104  and  were  stopped  when  Re^  reached  the  value  3.68  x  104,  which  is  the  value 

measured  by  Laderman  and  Demetriades  at  the  160-inch  station.  The  solutions 
show  that  even  over  a  streamwise  distance  of  about  805  inches  (corresponding  to 
the  increase  of  Re^  from  1,5  x  104  to  3.68  x  104  and  a  value  of  0  ■  0.29  in.  at 
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Re#  =  3.68  x  10^)  the  value  of  q  at  the  end  of  the  runs  is  still  dependent  upon 

the  initial  conditions  of  q  and2CH/q  .  Thus,  it  appears  likely  that  the  data 
of  Laderman  and  Demetriades  are  for  a  layer  relaxing  to  equilibrium  and  that 
because  the  throat  in  their  experiments  was  located  at  x  -  690  inches  (Figure 
22),  their  reported  value  of  'Cf  -  0.46  x  10“  3  is  probably  larger  than  the  equili¬ 
brium  flat  plate  at 'Re#  =  3.68  x  10^.  Because  the  flow  is  still  relaxing  and 

data  was  presented  at  only  one  streamwise  station  it  is  impossible  to  give  a 
precise  comparison  between  experiment  and  theory.  For  example,  it  is  i'3ssible 
to  match  the  experimental  values  of  0  =  0.29  in.  and  8*  =  2.29  in.  measured  at 
x  =  160  inches  with  various  pairs  of  values  of  q  and  2qj/q  at  the  intital 
station.  One  such  pair:  q  =  0.6  x  10“3  and(2qj/q)o  =  0.76  resulted  in  the 

values:  Cf  =  0.37  x  10“3,  o  =  0.29  in.,  8*  =  2.31  in.,  Se  a  5.53  in.  at 
Re #  a  3.68  x  10^,  while  the  values  measured  by  Laderman  and  Demetriades  were 
Cf  =  0.46  x  10-3,  0  =*.29  in.,  8*  =  2.29  in.,  and  8e  =  6.01  inch.  At  these  con¬ 
ditions  the  theory  predicted  that  wall  layer  extended  out  to  y  =  1.74  in.  where 
the  values  of  u/ue  and  H/He  were  computed  to  be  0.794  and  0.782,  respectively. 

(The  value  of  2  CH/q  at  Re^  B  3.68  x  10^  was  0.815).  The  profile  data  measured 
by  Laderman  and  Demetriades  gave  u/ue  =  0.77  and  H/He  =  0.74  at  y  *  1.74  in. 
However,  a  valid  comparison  would  require  data  to  be  obtained  at  two  and  prefer¬ 
ably  three  streamwise  stations  so  that  the  relaxation  process  could  be  duplicated. 
As  a  matter  of  reference,  for  theMe“  9.37  results  presented  in  Figure  22,  a 
perturbation  of  the  initial  Cf  by  20  percent  resulted  in  a  perturbation  of  cf  by 
5.5  percent  after  a  streamwise  relaxation  distance  of  160  inches  or  an  x/0o  m 
1.36  x  10-^ .  Consequently,  it  appears  that  virtually  all  data  obtained  on  side- 
walls^f  wind  tunnels  contain  upstream  relaxation  effects  ir»  Cf  as  well  as  in 
2  CH/Cf  and  some  of  these  data  may  be  far  out  of  fluid  dynamic  equilibrium. 

In  order  to  test  the  present  turbulent  boundary  layer  model  in  a  highly 
nonequilibrium  flow  situation,  a  solution  was  performed  for  flow  over  a  wavy 
wall.  Experiments  by  Bertram,  et.  al. ^  for  turbulent  boundary  layers  over  wavy 
walls  at  Mach  6  have  demonstrated  significant  surface  pressure  and  heating  vari¬ 
ations  (by  as  large  a  factor  as  2)  along  a  wave  train.  Maximum  heating  was 
measured  on  the  forward  face  of  each  wave  while  minimum  heating  was  obtained  on 
the  rearward  faces.  In  order  to  compute  the  turbulent  boundary  layer  develop¬ 
ment  over  waves  with  wave  length  of  the  order  of  the  boundary  layer  thickness, 
the  solution  was  obtained  as  a  quasi  boundary  layer-inviscid  flow  interaction. 

The  boundary  layer  at  the  start  of  the  wave  train  is  in  equilibrium  (dp/dx  =  0 ) 
for  Me  “  3.6  and  Re#  =  4.5  x  10-*.  The  equilibrium  boundary  layer  profile  at 

O 

these  conditions  was  used  as  the  initial  condition  for  a  rotational  characteris¬ 
tics  calculation  of  the  surface  pressure  distribution  along  the  waves.  The 
turbulent  boundary  layer  development  over  the  waves  was  then  computed  using  the 
surface  pressure  distribution  generated  by  the  characteristics  solution.  Figure 
23  shows  the  results  obtained  for  a  wave  train  with  wave  length  equal  to  the 
initial  boundary  layer  thickness  and  wave  amplitude  equal  to  0.004  times  the 
wave  length.  The  notation  max.,  min.  along  the  x  coordinate  refers  to  the  peaks 
and  troughs  of  the  surface  waves. 

The  response  of  the  momentum  thickness  to  the  periodic  surface  pressure 
variation  (or  edge  Mach  number  variation  shown  in  Figure  23)  is  characteristic 
of  supercritical  layers,  increasing  when  dp/dx  is  negative  and  decreasing  when 
dp/dx  is  positive,  giving  a  periodic  variation  of  0  with x  .  The  results  for  S’ 
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1  and  Se  are  also  periodic  with  roughly  the  same  phase  as  6.  As  expected,  the 

heat  flux  is  virtually  in  phase  with  the  surface  pressure  with  the  average  being 
;  somewhat  higher  than  for  a  flat  plate.  However,  the  average  skin  friction 

coefficient  is  less  than  the  flat  plate  value  and  the  streamwise  variation 
reveals ' local  minima  in  Cf  just  upstream  of  the  points  of  both  maximum  and 
;  minimum  in  surface  pressure.  This  peculiar  variation  of  C{  with  x  may  be  a 

result  of  the  supercritical  response  of  the  layer. 

.  I 

The  influence  of  the  additive  constant  B  in  the  generalized  compressible 
law  of  the  wall  on  the  turbulent  boundary  layer  over  a  sharp  20°  cone  at 
Mach  5.7  is  .-illustrated  in  Figure  24.  A  decrease  in  B ,  which  could  be  the  result 
,  of  an  increase  in  surface  roughness,-^  produces  an  increase  in  the  integral 

thicknesses  0,  8*  ,  and  Se  ,  and  an  increase  in  Cf  as  well.  Because  the  boundary 
layer  solutions '.shoyn  in  Figure  24  are  ^n  equilibrium  for  x~  10  inches,  with 
!  2CH/Cf  =  1,  the  Stanton  number  also  increases  with  decreasing  B.  For  rough 
walls 

—  _  I 

*  B,  f  B  -  4u]/u,  ’  .  | 

i  where  B  is  the  "constant"  for  smooth  walls.  For  roughness  densities  of  the  order 
4  to  10. 

i  '  1  1  ' 

A  uj  1  , 

i  •  — —  =  £„  Kyr/>'w)  , 

,ur  K  I  .  ! 

I  I 

*  ,  I  " 

j  where  yr  is  the 'roughness  height.  These  expressions  were  used  to  relate  yr  to 

'Bj  for  the  results  shown  in  Figure  24.  Evidently,  a  roughness  height  of  0.016 
•'  inch  produces  an  increase  in  by  a  factor  of  about  2.5  above  the  smooth  wall 

value. 

! 

!  i  ■ 

1  Calculations  showing  the  effect  of  a  sudden  expansion  on  turbulent  boundary 
layer  cheating  are  compared  with  some  measurements  in  Figures  25,  26,  and  27. 

I  The  experiments  were  conducted  at  a  free  stream  Mach  number  of  5  in  a  hypersonic 

wind  tunnel  at  the  U.S.  Natal  Ordnance  Laboratory.  The  measurements  were  per¬ 
formed  on  a  sharp  50-8  degree  biconic  model  and: are  reported  in  Reference  27. 
j  The  results  presented  in  Figures  25  through  27  are  for  three  values  of  the  free 
stream  unit  Reynolds  number.  As  indicated  by  the  heat  flux  levels  reported  in 
l  Reference  27,  the  boundary  layer  was  fully  turbulent  over  at  least  half  of  ,the 

50  degree  cpne  for  all  three  Reynolds  numbers. 

» 

.  At  the  50  and  8  degree  cone  junction  both  experiment  and  theory  show  a  sud¬ 

den  decrease  in  heating,  however,  the  data  shows  that  for  a  certain  distance 
! after  the  expansion  the  heaping  falls  well  below  theoretical  predictions.  At 
the  two  highest  values  of  -Reynolds  number  (Figures  25  and  26)  the  heating  drops 
to  about  30  percent  of  the  theoretical  prediction  at  the  end  of  the  expansion 
but  after  a  certain  distance  downstream  (about  65  times  the  boundary  layer  thick- 
!  ness  'before  the  expansion  for  the  case'  shown  in  Figure  26)  it  increases  to  about 

!  'the  level  predicted.  At  the  lowest  Reynolds  number  (Figure  27)  the  heating 

I  *  . 

*A  continuous  turbulent  boundary  Iftyer  sdlution  over  the  sudden  expansion  was  obtained  t>,  assuming  that  the  complete 
expansion  from  0C  =  50° (Me  =  1.15)  to  <?c  =  8°(MC  =  2.73)  takes  place  over  a  streamwise  distance  of  0.16  in.  or  about 
•  5  local  boundary  layer  thicknesses  (8  =  0.032  in.  ahead  of  the  expansion). 

1  *  I  i  C 
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Figure  26  RELAM1NARIZING  FLOW:  HEAT  TRANSFER  Figure  27  RELAMINARIZING  FLOW:  HEAT  TRANSFER 

DISTRIBUTION  ON  A  SHARP  BICONIC  DISTRIBUTION  ON  A  SHARP  BICONIC 


drops  suddenly  at  the  expansion,  as  in  the  other  two  cases,  but  then  continues 
to  decrease  at  a  fairly  rapid  rate  along  the  8  degree  cone.  At  the  last  data 
point  (x7R  =*  1.93)  the  heating  is  only  about  25  percent  of  that  predicted  by  the 
theory.  Apparently,  the  decrease  in  heating  across  the  expansion  is  controlled 
by  the  pressure  gradient  for  these  cases  because  q /q^  at  the  start  of  the  8 
degree  cone  is  approximately  the  same  for  all  three  Reynolds  numbers. 

The  behavior  of  the  heating  data  in  these  experiments  suggests  that  the 
turbulent  boundary  layer  on  the  50  degree  cone  partially  or  even  completely 
relaminarized  in  the  sudden  expansion  at  the  junction  of  the  50  and  8  degree 
cones.  The  abrupt  increase  in  q  at  some  distance  downstream  of  the  expansion  at 
the  two  highest  Reynolds  numbers  suggests  a  retransition  of  the  laminarized 
layer.  Evidently,  at  the  lowest  Reynolds  number  retransition  occurred  downstream 
of  the  last  heating  gage  so  that  the  flow  remained  laminarized  for  at  least  90 
boundary  layer  thicknesses  (based  on  the  layer  thickness  just  upstream  of  the 
expansion) . 

•  Several  experiments  28,29  on  turbulent  boundary  layer  laminar ization  in  in¬ 
compressible  flows  have  demonstrated  that  during  a  strong  expansion  the  sadden 
reduction  of  turbulence  production  in  the  wall  layer  produces  a  departure  of  the 
velocity  profiles  from  the  law  of  the  wall.  When'u+  is  plotted  against  y'*  the 
profiles  shift  progressively  upward  away  from  the  law  of  the  wall  an  the  flow 
proceeds  through  a  sufficiently  strong  expansion.  Launder  and  Jones^®  have  sug¬ 
gested  t’.it  this  upward  shift  in  the  velocity  profile  and  the  simultaneous  in¬ 
crease  in  the  laminar  sublayer  thickness  can  be  characterized  by  an  increase  in 
the  parameter  a  in  the  relation  for  the  variation  of  mixing  length  through  the 
sublayer  and  transition  layer,  Eq.  (55).  As  shown  in  Figure  3,  an  increase  in  a 
produces  an  increase  in  the  additive  "constant"  B  in  the  law  of  the  wall  so  that 
one  of  the  observed  features  of  laminarization  can  be  simulated  by  permitting  B 
to  increase_with  7  through  an  expansion.  Figures  26  and  27  show  results  obtained 
by  letting  B  increase  linearly  with  7  through  the  expansion  (a7  =  0JL62  in.) 
from  the  initial  value  of  6.1  (for  Hw/He  =  0.4)  to  selected  values  (B2  “  18,  25 
and  3i)  at  the  end  of  the  expansion.  Along  the  8  degree  cone  following  the 
expansion,  B  a  B2  a  constant.  These  results  demonstrate  that  a  decrease  in  heat¬ 
ing  below  the  turbulent  prediction  is  consistent  with  an  increasing  B  (or  _ 
increasing  a)  in  the  expansion.  By  modeling  the  variation  of  a  rather  than  B 
through  strong  expansions,  and  by  including  the  contribution  of  the  laminar 
sublayer  in  the  relations  for  mass,  momentum  and  energy  flux  in  the  wall  layer, 
it  should  be  possible  to  reproduce  even  more  closely  some  of  the  features  of 
laminarization  and  retransition  displayed  by  the  data  in  Figures  25  through  27. 

7.2  Results  for  Uniform  and  Discontinuous  Distributions  of  Surface  Mass  Injection 

The  predicted  variation  of  Cf  with  Re^  for  various  blowing  rates  in  incom^ 
pressible  flow  is  shown  in  Figure  28.  The  results  sho^that  at  a  given  Re^  ,  Cf 
decreases  steadily  to  zero  as  M  increases.  Values  of  Cf  as  low  as  10**15  have 
been  computed  for  M  =  0.05  and  Reg  ^  10$  so  that  with  the  two-layer  model  there 
is  ho  evidence  of  an  Emmons  type  singularity  for  massive  blowing. 

Although  the  experimental  data  for  C,  with  surface  mass  injection  is  highly 
scattered,  Squire®  has  attempted  a  best  fit  of  data  obtained  by  McQuaid30, 
Jerondr^1  and  Squire®  at  Me  -  0,  2.5  and  3.5,  which  is  compared  with  the  theory 
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in  Figure  29.  Following  Squire  and  others,  the  Cf/Cf0  ratio  is  plotted  versus 
the  usual  injection  parameter*  2M /C^  .  Both  the  experiments  and  theory  indicate 
that  the  decrease  in  Cf  at  a  given  value  of  the  blowing  parameter  is  less  severe 
as  the  Mach  increases. 

One  possible  reason  for  the  large  scatter  in  Cf  data  with  blowing  is  shown 
in  Figure  30,  where  we  have  compared  predictions  with  the  data  of  Danberg9  at 
Mach  6.7.  When  plotted  in  the  standard  way,  i.e.,  Cf/Cfo  versus  2M/Cf0  ,  the 
theory  reveals  a  significant  variation  with  the  local  Re# ,  which  apparently  has 
not  been  accounted  for  in  past  presentations  of  data.  The  theory  also  predicts 
a  variation  in  the  Sf/Cf  ratio  with  wall  to  free-stream  stagnation  enthalpy 
ratio,  as  shown  in  Figure  30. 

_  Thg  effect  of  surface  mass  injection  on  the  form  parameter  H  = 8*/0  ,  where 
8  and  0  are  the  displacement  and  momentum  thicknesses  in  transformed,  Howarth  co¬ 
ordinates  is  shown  in  Figure  31.  The  data  extends  over  a  Mach  number  range  of  0 
to  6.2,  and  was  compiled  by  Fernandez  and  Zukoski32.  Results  of  the  two-layer 
model  were  obtained  for  Me  =  0  and  are  seen  to  agree  quite  well  with  the  data  up 
to  the  largest  value  of  M  (0.013)  for  which  there  are  measurements  at  Me  <=  0. 

For  M  >  0.015  the  theoretical  curve  diverges  rapidly  from  the  data,  which,  for 
these  large  blowing  rates  was  all  obtained  for  Me  =  2.5.  Thus,  it  is  not  clear 
whether  the  divergence  of  the  curve  from  the  data  signals  the  failure  of  the  two- 
layer  model  for  large  blowing  rates  or  a  breakdown  in  the  Howarth  transformation 
in  collapsing  all  compressible  profiles  to  a  constant  density  form. 

An  increase  in  the  physical  8*/0  ratio  with  increasing  blowing  is  predicted 
at  all  Mach  numbers.  For  example,  Figure  32  shows  a  comparison  of  results  at 
Mach  6.7  with  data  obtained  by  Danberg9,  while  Figures  33  and  34  show  comparisons 
with  the  streamwise  variations  of  0  and  8 *  measured  by  Jeromin31  at  Mach  2.5  and 
3.5,  respectively. 

Figure  35  gives  a  comparison  of  predicted  momentum  thickness  variations  with 
the  measurements  of  Martellucci33  for  uniform  injection  through  a  7.25-degree 
half-angle  cone  at  =  8.  All  solutions  were  started  at'x*5  8.4  inches,  which 
was  the  station  at  which  injection  began  on  the  model.  This  station  also  corres¬ 
ponded  approximately  to  the  point  where  boundary  layer  transition  began  in  the 
experiments,  having  evidently  been  tripped  by  the  injection,  while  the  end  of 
transition  was  located  at  about  x  «  16.5  inches.  Two  solutions  were  performed 
for  M  =  0,  starting  with  =  0.003  and  0.004  inch,  to  test  whether  the  down¬ 
stream  0  variations  would  be  strongly  influenced  by  the  initial  conditions  for 
Cf,  and  8e  ,  which  were  not  measured  at  x'  =  8.4  inches  and,  therefore,  had  to 
be  assumed.  The,  two  solutions  for  M  =  0  are  seen  to  converge  rapidly  to  a  more 
or  less  unique  0  variation  with 7 ,  indicating  that  the  solutions  have  relaxed  to 
equilibrium  for  'x  >  25  inches.  It  is  assumed  that  the  solutions  for  M  £  0  have 
also  achieved  equilibrium  and  are  virtually  independent  of  the  starting  conditions 
for  x*  >  25  inches. 

Figure  35  shows  that  for  M  =  0  and  the  two  lowest  injection  rates  the  theory 
predicts  values  of  0  well  below  the  data,  while  for  the  two  highest  injection 
rates  the  results  are  in  good  agreement.  Since  the  measurements  were  all  obtained 


*  Examination  of  the  equations  for  the  two-layer  model  (Appendix  1)  shows  that  the  proper  injection  parameter  is  2M/Cf . 
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Figure  33  VARIATIONS  OF  6  AND  5*  AT  Me  =  2.5  Figure  34  VARIATIONS  OF  8  AND  6*  AT  Me  =  3.5 


with  the  porous  sintered  steel 
model,  the  wide  disagreement  be¬ 
tween  the  experiments  and  theory 
for  M  *  0  is  apparently  caused 
by  surface  roughness  (see,  for 
example.  Figure  24).  But  as  M 
increases,  producing  larger 
boundary  layer  thicknesses  com¬ 
pared  to  the  roughness  height, 
the  influence  of  surface  rough¬ 
ness  evidently  disappears.  These 
results  demonstrate  that  the 
combined  effects  of  mass  injection 
and  surface  roughness  can  be 
important  and  will  require 
further  study  to  produce  reliable 
boundary  layer  predictions, 
especially  for  very  cold  walls 
where  effects  of  surface  rough¬ 
ness  are  magnified  and  where 
there  is  little  or  no  data  on  the 
effect  of  Hw/He  and  M  on  the  addi¬ 
tive  constant  b  in  the  law  of  the 
wall.  It  is  also  quite  likely 
that  roughness  effects  are  re¬ 
sponsible  for  must  of  the  scatter 
in  Cf  and  Cjj  data  with  surface 
mass  injection. 


,  „  ^  ,  Figure  35  MOMENTUM  THICKNESS  DISTRIBUTIONS  ON  A 

Solutions  for  discontinuous  SHARP  CONE  AT  MACH  8 

blowing  distributions  in  in¬ 
compressible  flows  are  compared  with  some  measurements  by  McQuaid^O  and  Simpson^ 
in  Figures  36  and  37.  For  the  step  decrease  in  blowing  data  on  McQuaid,  the  solu¬ 
tion  was  performed  by  constructing  a  continuous  ramp  function  for  M,  which  de¬ 
creased  from  3.4  x  10“3  to  zero  over  a  streamwise  distance  equal  to  the  local 
momentum  thickness.  The  initial  abrupt  increase  in  Cf  shows  that  at  the  first 
station  the  theoretical  solution ^is  out  of  equilibrium,  with  Cf  being  below  the 
theoretical  equilibrium  value.  Cf  then  reaches  the  equilbrium  value  and  then  de¬ 
creases  with  x  to  the  point  where  the  injection  stops. 


A  continuous  ramp  function  was  used  to  increase  M  from  zero  to  0.002  over 
a  streamwise  distance  equal  to  the  local  0  in  the  solution  given  in  Figure  37. 

The  solid  curves  show  the  behavior  of  the  solution  for  M  *  0.002  downstream  of 
x  a  37.5  inches,  while  the  dashed  curves  show  the  behavior  if  the  blowing  had 
remained  equal  to  zero. 

Velocity  profiles  with  surface  mass  injection  (see  Section  6.0)  are  compared 
with  incompressible  data  of  McQuaid^O  in  Figure  38  and  with  data  obtained  by 
Jeromin^l  at  Me  =  3.5  in  Figure  39.  In  both  cases  a  ,  which  is  determined  by 
matching  B  computed  from  Eq.  (59)  with  the  empirical  relations  for  B  given  in 
Section  3.0,  decreases  steadily  with  increasing  M  (see  Figure  1).  At  some 
sufficiently  large  value  of  M  ,  depending  on  Hw/He  and  Mc ,  a  ■*  0,  which  implies 
that  the  laminar  sublayer  vanishes.  For  incompressible  flow  a  vanishes  at 
M  =»  0.01  so  that  in  order  to  maintain  a  -  0  for  M  =>  0.01  B  would  have  to  increase 
with  increasing M.  This  conjectured  behavior  of  B  for  large  blowing  is  just 
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DISCONTINUOUS  DECREASE  IN  BLOWING,  Figure  37  DISCONTINUOUS  INCREASE  IN  BLOWING, 

INCOMPRESSIBLE  FLOW  INCOMPRESSIBLE  FLOW 


Figure  38  VELOCITY  PROFILES  WITH  INJECTION,  Me  =  0  Figure  39  VELOCITY  PROFILES  WITH  INJECTION,  M( 


opposite  to  the  variation  of  B  with  M  which  was  established  from  the  data  of 
Squire  and  Danberg  for  small  blowing  at  supersonic  speeds. 

As  M  increases  above  0.01,  the  velocity  profiles  obtained  using  the  law  of 
the  wall  in  the  inner  layer,  rather  than  Eqs.  (57)  and  (58),  become  highly  in¬ 
flected,  with  most  of  velocity  increase  taking  place  in  the  outer  part  of  the 
boundary  layer.  This  behavior  was  also  found  in  the  analyses  of  turbulent 
boundary  layers  with  massive  blowing  in  References  35  and  36  and  in  the  experi¬ 
ments  reported  in  Reference  32. 

7.2.1  Ablation  and  Transpiration  Cooling:  Solutions  for  Flows  with  the  Mass 
Injection  Distribution  Determined  by  an  Energy  Balance  at  the  Surface 

A  crucial  feature  of  flows  with  ablation  or  transpiration  cooling  is  the 
interaction  mechanism  between  the  rate  at  which  gas  is  sublimated  or  evaporated 
into  the  boundary  layer  and  the  heat  flux  which  reaches  the  solid  or  liquid 
surface.  In  this  situation  the  mass  injection  rate  is  determined  by  a  local 
energy  balance  at  the  surface  and  the  blowing  parameter  M  becomes  an  additional 
dependent  variable  in  the  system  of  simultaneous  differential  equations. 

If  Hj  is  the  effective  heat  of  ablation  or  evaporation,  which  includes  the 
heat  capacity,  i.e.,  H;  «  Hv  +  Cp  AT  ,  where  Hv  is  heat  of  sublimation  or 
evaporation,  then 

=  (65> 


For  negligible  gas  radiation  and  negligible  conduction  in  the  liquid  or  solid, 
qw  is  the  convective  heat  flux.  Combining  Eqs.  (37)  and  (65)  gives 

M  =  He/Hj  (1  -  Hw/He)  CH  (66> 

From  Eqs.  (25)  and  (37) 


CH  = 


Hw/He 

-  BCt 

(1  -  Hw/He)  1 


(67) 


Substituting  Eq.  (67)  into  (66)  and  differentiating  with  respect  to  the  stream- 
wise  distance  x  gives  (assuming  H;  is  independent  of  x ) 


dM  Hw  ["  dB  dCf‘ 
dx”  =  Hj”  f  "dsT  +  B  ~dx"_ 


(68) 


By  using  Eq.  (66)  as  an  additional  initial  condition,  and  by  adding  Eq.  (68)  to 
the  system  of  simultaneous  differential  equations,  the  coupled  "feedback"  problem 
of  the  effect  of  blowing  on  the  heat  flux  can  be  computed  as  part  of  the  down¬ 
stream  integration  scheme.  The  use  of  Eq.  (68)  obviates  having  to  perform 
lengthy  iterations  for  the  proper  matched  boundary  conditions  on  M  and  at 
each  point  along  the  surface  so  that  the  computing  time  is  about  the  same  as  if 
M  is  a  specif! function  of  ~  . 
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Figure  40  shows  some  results  computed  with  the  coupled  energy  balance  at  the 
surface  for  a  0.25-inch  radius  spherical  nose  at  hypersonic  flow  conditions.  The 
stagnation  pressure  behind  the  normal  shock  is  40  atm,  the  free-stream  stagnation 
enthalpy  is  2,600  Btu/lb  and  the  wall  temperature  is  1000°  R,  resulting  in 
Hw/He  m  0.092.  At  these  conditions,  y  *  1.23  and  o  —  0.5  for  equilibrium  air. 

The  Mach  number  distribution  around  the  sphere  was  computed  assuming  isentropic 
flow  and  a  static  pressure  distribution  given  by 

p/Pt2  »  0.77  cos2d  +  0.23 

The  solutions  were  all  started  at  10  degrees  from  the  stagnation  point  using 
more  or  less  arbitrary  initial  conditions  for  “o  ,  %  ,  Cf  ,  CH .  Solutions  ob¬ 
tained  using  different  sets  of  initial  conditions  gave  results  for  Cf  and  Cjj 
that  were  virtually  independent  of  the  initial  conditions  for  angles  greater 
than  25  degrees  from  the  stagnation  point.  Figure  40  shows  the  effect  of  the 
parameter  He/Hi  on  distributions  of  the  blowing  parameter  M  and  the  nondimensional 
heat  flux.  As  the  parameter  He/H;  increases  the  rate  at  which  material  is  ablated 
or  evaporated  into  the  boundary  layer  increases,  resulting  in  a  net  decrease  in 
the  heat  flux  which  reaches  the  gas-solid  or  gas-liquid  interface. 


82-3404  */Rn 


Figure  40  BLOWING  AND  HEATING  DISTRIBUTIONS  AROUND 
A  SPHERICAL  NOSE  WITH  COUPLED  ENERGY 
BALANCE  AT  THE  SURFACE 
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8.0  CONCLUDING  REMARKS 


i 
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All  results  reported  In  this  report  have  been  obtained  tiy  evaluating  the 
turbulent  production  integral  over  the  outer  layer  using  a  cubic  variation  of  c 
between  ym  and  Se  (w=  3  in  Eq.  (46)).  Stable  solutions,  in  the  sense  that  the 
power-law  exponent  n  in  the  outer  layer  velocity  profile  remains  bounded,  have 
been  obtained  for  negative  and  positive  pressure  gradients,  with  cold  and  ■ 

adiabatic  walls,  up  to  Mach  10.  It  is  really  quite  remarkable  that  the  evalua¬ 
tion  of  the  production  integral  with  a  simple  cubic  variation  of  (  produces 
stable  solutions  for  such  \  wide  range  of  flow,  conditions.  For  incompressible 
flow  in  zero  pressure  gradient,  for  example,  a  constant  value  of  the  production 
integral  aQ  equal  to  0.18  gives  a  downstream  solution  of'  the  type  n  -*  o°  ,  while  a 
value  a0  =  0.09  gives  a  solution  in  which  n  -  0 .  Moreover,  these  'instabilities 
occur  after  a  streamwise  distance  of  only  about  10  boundary  layer  thicknesses. 
Despite  this  rather  sensitive  dependence  of  the  streamwise  variation  of  n  on  the 
magnitude  of  aQ  ,  the  values  of  Cf  along  the  integral  curves  differ  by  less  than  • 
one  percent  for  aQ  =  0.18  and  a0  =  0.09.  Clearly,  then,  the  stability  of  the 
interaction  between  the  inner  and  outer  layers  is  much  more  sensitive  to  the 
level  of  turbulence  production  in  the  outer  layer  than  is  the  streamwise  develop¬ 
ment  of  the  layer  as  a  whole.  In  other  words,  as  long  as  a  «a\  of  calculating 
aQ  can  be  found  which  keeps  the  inner-outer  layer  interaction  stable  ( n  finite) , 
the  results  for  &  (~),  'o  Cx) ,  etc.  should  be  within  a  few  percent  of  the  results 
of  all  other  stable  solutions  obtained  with  slightly  different  methods  of  cal¬ 
culating  a0  .  This  is  true  whether  a0  is  computed  using  eddy  viscosity  or  other, 
more  sophisticated,  stress  models. 

Solutions  for  boundary  layers  in  strong  expansions  have  exhibited  instabili¬ 
ties  of  the  type  n  <*>  for  w  =  3  so  that  the  cubic  variation  of  c/rm  obviously 
cannot  be  considered  "universal."  Because  for  this  type  of  instability  (n -  °°  ) 
aQ,  and  hence  o,  must  be  decreased  to  produce  stable  solutions,  a  quadratic  or 
even  linear  variation  of  f/fm  across  the  outer  layer  may  be  required.  An 
example  of  the  unstable  behavior  in  a  strong  negative  pressure  gradient  is  shown 
in  Figure  41.  The  calculations  were  performed  for  a  flow  in  which  the  inviscid 
Mach  number  increases  linearly  from  1.5  to  4.5  over  a  streamwise  distance  equal 
to  660  times  the  initial  boundary  layer  thickness.  Solutions  were  performed  for 
(d«  1,  2  and  3,  starting  with  an  initial  boundary  layer  in  equilibrium  for  zero 
pressure  gradient  at  Re^  =  10^.  As  indicated  in  the  figure,  n  grows  to  10^  for 

(o  =  3  and  2  at  x/S  =  390  and  430,  respectively.  Both  solutions  are  unstable 
co  ' 

with  n  increasing  exponentially  at  these  values  of  x  .  The  solution  for  n  =  1  is 

stable  throughout  the  expansion,  with  n  a  27.6  at  the  end  of  the  run.  ’Despite 
the  wide  difference  in  the  stability  characteristics  of  the  solutions,  Figure  41 
shows  that  the  maximum  difference  in  Cf  along  the  integral  curves  for  the  three 
values  of  «  is  only  about  4  percent.  A  stable  solution  is  founi  by  lowering  the 
value  of  w  (decreasing  a0 )  from  3  to  1  but  the  effect  of  this  operation  on  the  I 
boundary  layer  solution  is  practically  negligible.  Consequently,  it  appears 
that  the  two-layer  model  can  provide  quantitative  results  for  the  streamwise 
variation  of  the  turbulence  production  in  the  outer  layer  simply  on  the  basis  of 
the  stability  of  the  inner-outer  layer  interaction.  Solutions  obtained  using 
finite  difference  methods11*12  and  prescribed  lateral  variations  of  mixing 
length  or  eddy  viscosity  have  also  shown  that  the  overall  boundary  layer  develop¬ 
ment  is  rather  insensitive  to  assumed  levels  of  intermittency  at  the  outer  edge. 
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However,  these  methods  are  incapable  of  providing ’much  insight  into  the  inner 
!  outer. layer  interaction,  at  least  as  theytare  presently  formulated. 


v  Figure  41  SKIN  FRICTION  VARIATIONS  IN  A  NEGATIVE 
PRESSURE  GRADIENT  WITH  w=  1.  2.  AND  3 
IN  P IRBULENT  PRODUCTION  INTEGRAL,  aQ 
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APPENDIX  1 


A 


MATCHING  CONDITIONS 


M 


In  addition  to  the  law  of  the  wall  for  Um  and  the  Crocco  integral  for 
H/H^,  matching  relations  for  (du/dy  )m  and  (dH/dy )  m  are  required  to  determine 
the  exponents  n  and  m  of  the  power-law  profiles  in  the  outer  layer.  These  are 
given  below  along  with  the  expressions  for  rm ,  (q  +  ur  )m  and  m  the  mass  flux 
in  the  inner  layer.  These  expressions  are  applicable  to  a  boundary  layer 
growing  into  a  rotational  (variable  entropy  but  constant  stagnation  temperature) 
inviscid  flow. 


From  continuity  of  (du/dy  )  at  ym: 

y  + 1 

Po  2  ( y -  1 )  3/2  Cf  J/2  6  l'2 

n  =  ——  ( 1  +  «n„ )  •  - 
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From  continuity  of  (dH/dy)  at  ym: 
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where 
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Stress  at  ym  (obtained  by  differentiating  J2  and  inserting  results  into 

relation  for  r  : 
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Heat  flux  at  ym  (from  equation  for  (q  +  ur  )m  ): 


-62- 


where 


'll  = 


*  cf1/2 


f  yffl 


Po  .L  <hm/hw>1/2 

1-3  y 

+  (1  +me)  2  (y-  D  Jj  B 


01  Cf1/2 


1  - 


Um  (B-2A2Um) 
2<hoAw> 


20O2  C* 


+  h  ym  M 


*12 


Ps 

"13  n  ~ 


P°  L  (hm/hw >  1/2 

$1  I'm 


+  02  0O2  Cf 


P°  L  2(hm/hw)1/2Cf1/2 
l-3y 


l  - 


0O  Cf 


*14 


“15 


“16  = 


Po  L 


fl,  Cr1'^  U  2  V 

(1  +  mc)2(y_1)  U1um  +  J2)  +  — f  ■  Xm  . 

2(hm/hw)3/2  J 


+  020o  ym  -  0 2  ymMum/c; 

1/2  77  2 


Ps 


P  2  Vm  ~  ' 


0i  Um2  y„ 


P2  \  'm  0o‘ 


M. 


1  - 


/VA  me 

^  y- 1  J  ( 1  +  rn. ) 


01  Cfl/2  Um3  ym  (He/Hw)(y-l)Me 
2(hm/hw)3/2  (l+me)2 


He) 1/2  BCf1/2  Um  ym  ft 


y 

y-  1 


*17 


KMe  0^/1^)  1/2 

+  Pz  cf  Po2  y, 


( 1  +  me  ) 


01  Cf1/2  Um  y® 


I1-! 

^  2y  \  me 

r 

y  y-  ly  (1  +  me  )  1 

_ 

' 

<hm/hw)1/2 


1  ~3y 


“18  D  2(ps/p0)®e  (l+me)  2^y-D 

-(y  +  D 

“19  "  (Ps/Po)  (1  +me)  2(y-l)  (1^/Hg)  B 

01  -  («w/He)1/2  ( 1  +  m# )  y~  1  (Bft0/K) 

-<y  + 1) 

02  -  2  (I^/Hg)  (1  +me)  2<y-D  (B/K2) 


-63- 


The  mass  flux  in  the  inner  layer  is 
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Also,  from  the  definition  of  B 
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The  velocity  ratio  Um  along  the  match  point  y  =  ym  is  found  fron  the  law  of 
the  wall,  Eq.  (26).  Differentiating  this  equation  with  respect  to  x  results  in 
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and  5(0, k)  ±s  the  elliptic  integral  of  the  second  kind 
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Differential  equations  for  dHm/dx,  dn/dx,  dm/dx  and  (l/p0  Ue)  (<Wdx)  are  also 
required  as  part  of  the  system  of  equations  to  be  solved,  and  are  obtained  by 
differentiating  Eqs.  (A  1-1),  (A  1-2),  (A  1-3),  and  (A  1-7),  respectively. 

Since  all  the  differential  equations  obtained  from  the  matching  conditions 
are  somewhat  lengthy,  a  system_of  jredundant  checks  has  been  built  in  to  the  com¬ 
puter  program.  The  values  of  Um ,  Hm,  n,  m  and  m/p0  Ue  calculated  from  the  system  of 
simultaneous  differential  equations  are  compared  at  each  station  wirh  the  values 
obtained  from  the  associated  algebraic  equations  for  these  quantities. 

Two  additional  redundant  checks  are  provided  by  overall  momentum  and  energy 
integrals  over  the  whole  boundary  layer. 


The  momentum  integral  equation  is 
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and  the  energy  integral  equation  is 
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where 


which,  when  evaluated  using  the  two-layer  model,  is 
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Values  of  9  and  obtained  by  integrating  Eqs.  (A  1-10)  and  (A  1-11)  are 
compared  with  values  computed  using  Eqs.  (51)  and  (A  1-12)  at  each  station. 
Special  computer  runs  have  been  made  using  double  precision  computations  with 
exceptionally  tight  error  control  in  the  simultaneous  linear  equation  solver  and 
in  the  Runge-Kutta  integration.  For  these  runs  the  algebraic  relations  for  vari¬ 
ables  along  the  match  point  and  for  9  and  agreed  with  the  results  computed 
from  the  differential  equations  to  within  nine  significant  figures. 
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Integral  functions  for  integral  moment  equations  in  outer  layer 
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APPENDIX  3  1 

ENTRAINMENT  AND  OTHER  SWALLOWING  RELATIONSHIPS 

l 

For  flow  around  blunt  bodies  at  angle  of  attack  the  streamline  locations  and 
metrics  are  calculated  using  the  method  developed  by  Fannelop  and  Waldman^.  In 
this  calculation  the  local  entropy  (or  stagnation  density  ratio  p's/p0  )  is  related 
to  the  mass  flow  entrained  by  the  boundary  layer  between  neighboring  inviscid  . 
streamlines.  If  me  is  the  mass  flow  entrained  by  the  boundary  layer  between  two 
inviscid  streamlines,  which  are  separated  by  a  distance  r  , 


^  =  tf  p'udy’  -  J  pw‘vwrd'? 
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(A  3-1) 


It  is  assumed  here  that  no  mass  is  injected  or  ablated  into  the  boundary 
layer  upstream  of  f  »  0  which  is  the  station  where  the  turbulent  boundary  layer 
begins.  Eq.  (A  3-1)  can  be  written 
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The  local  Reynolds  number  base  on  momentum  thickness  is 
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while  the  local  heat  flux  and  shear  stress  along  the  surface  are 
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